Crystal Growth and Physical Properties of SrCu 2 As 2 , SrCu 2 Sb 2 and BaCu 2 Sb 2 
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We report the growth of single crystals of SrCu2As2, SrCu2Sh)2, SrCu2(Aso.84Sbo.i6)2 and 
BaCu2Sb>2 using the self- flux technique and their structural, magnetic, thermal and transport prop- 
erties that were investigated by powder x-ray diffraction (XRD), magnetic susceptibility x> specific 
heat C p and electrical resistivity p measurements versus temperature T from 1.8 to 350 K. Ri- 
etveld refinements of XRD patterns for crushed crystals confirm that SrCu2As2 crystallizes in the 
ThCr2Si2-type body-centered tetragonal structure (space group IA/mmm) and SrCu2Sb2 crystal- 
lizes in the CaBe2Ge2-type primitive tetragonal structure (space group PA/rnnm). However, as 
reported previously, BaCu2Sb2 is found to have a large unit cell consisting of three blocks. Here 
a ThCr2Si2-type block is sandwiched between two CaBe2Ge2-type blocks along the c-axis with an 
overall symmetry of I4/mmm, as reported, but likely with a monoclinic distortion. The \ data of 
all these compounds are diamagnetic and reveal nearly T-independent anisotropic behavior. The \ 
of SrCu2As2 is found to be larger in the afc-plane than along the c-axis, as also previously reported 
for pure and doped BaFe2As2, whereas the \ values of SrCu2Sb2 and BaCu2Sb2 are larger along 
the c-axis. This difference in anisotropy appears to arise from the differences between the crystal 
structures. The finite values of the Sommerfeld linear specific heat coefficients 7 and the T depen- 
dences of p reveal metallic character of all four compounds. The electronic and magnetic properties 
indicate that these compounds are sp metals with Cu in the nonmagnetic 3d 10 electronic configura- 
tion corresponding to the oxidation state Cu +1 , as previously predicted theoretically for SrCu2As2 
by D. J. Singh [Phys. Rev. B 79, 153102 (2009)]. We present a brief review of theoretical and 
experimental work on the doping character of transition metals for Fe in BaFe2As2. The As-As co- 
valent interlayer bond distances in the collapsed-tetragonal (Ca,Sr,Ba)Cu2 As2 compounds are much 
shorter than the nonbonding As-As distances in BaFe2As2- Thus the electronic character of the Cu 
and the strength of the As-As interlayer bonding are both expected to drastically change between 
weakly Cu-substituted BaFe2As2 and pure BaCu2As2, perhaps via a first-order lattice instability 
such as a miscibility gap in the Ba(Fei_ a ,Cu :E )2As2 system. 

PACS numbers: 74.70.Xa, 72.15.Eb, 65.40.Ba, 74.70.Dd 



I. INTRODUCTION 



The observation of superconductivity with transition 
temperatures T c up to 38 K in 122-type iron pnic- 
tides (e.g., J 4i_. r K; E Fe2As2 compounds, A = Ca, Sr, 
Ba, and Eu) 1-5 having the layered ThCr2Si2-type struc- 
ture, shortly after the discovery of superconductivity 
with T c up to to 55 K in 1 Ill-type oxypnictides (e.g., 
LnFeAsOi-ajFj; compounds, Ln = La, Ce, Nd and 
Sm), 6 ~ 10 triggered research activities worldwide and it 
is now believed that the iron pnictide compounds rep- 
resent a new class of high-T c superconductors. 11-13 Su- 
perconductivity in these iron arsenides lies in close prox- 
imity to antiferromagnetic (AFM) itinerant spin density 
wave (SDW) transitions that can be easily suppressed 
by substitutions at the A, Fe and/or As sites. 11 The 
generic phase diagram for the emergence of superconduc- 
tivity upon suppression of the SDW transition in iron 
arsenides is qualitatively similar to that of the high-T c 
cuprates, although with the important difference that 
the cuprate parent compounds are local moment antifer- 
romagnetic insulators whereas the iron arsenide parent 
compounds are itinerant SDW semimetals. 11-17 Interest- 
ingly, as noted above superconductivity in iron arsenides 
can be induced by substitutions on the Fe site by other 
transition metals, which is quite different from the high- 



Tc cuprate superconductors where such substitutions at 
the Cu site never induce superconductivity and indeed 
can lead to a rapid suppression of the superconductivity. 

In the fiducial 122-type compound BaFe2As2, the for- 
mal oxidation states of the atoms are assigned as Ba +2 , 
Fc +2 and As -3 so the Fe atoms are formally in the 3c? 6 
electronic configuration as also occurs in the 1111-type 
parent compounds. 11 It has been observed that partial 
Co and Ni substitutions at the Fe site in BaFe2As2 in- 
duce superconductivity 18-20 with T c up to 25 K whereas 
no superconductivity is induced by Mn (Refs. 21-23) 
or Cr (Refs. 24, 25) substitutions. These findings sug- 
gest that the superconductivity might have some rela- 
tion with the average number of 3c? conduction elec- 
trons of the transition metal atom. That is, hole-doping 
with a smaller number of 3c? electrons (Mn +2 , Cr +2 ) 
than that of Fe +2 does not induce superconductivity, 
whereas electron-doping with more 3c? electrons than 
Fe +2 (Co +2 and Ni +2 ) does induce superconductivity. 
From this point of view divalent copper Cu +2 with the 
3c? 9 electronic configuration and three more c?-electrons 
than Fe should be a strong electron dopant for iron ar- 
senide superconductors. However, Cu-doping for Fe in 
BaFe2As2 has been found to yield strongly suppressed 
T c < 2 K in Ba(Fei_ a: Cu a: )2As2, and then only in a 
very limited concentration range near x = 0.044, even 
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though the Cu doping suppresses the structural/ SDW 
transition of the parent compound BaFe 2 As2. 26 Simul- 
taneous Co- and Cu-doping for Fe in BaFe2As2 leads 
to superconductivity with higher values of T c for y > 
in Ba(Fei_ a; _j ) Coj / Cu 2: )2As2 compounds than for y = 
(Ref. 26). On the other hand, the application of pressure 
to AFe2As2 (A — Sr or Ba) compounds or isoelectronic 
substitutions of Ru for Fe or P for As in BaFe2As2, which 
nominally do not result in charge doping, can also induce 
high-T c superconductivity. 11 

A +2 oxidation state of copper also plays a key role 
in high-T c cuprates such as in the antiferromagnetic in- 
sulator parent compound La2Cu04 with a Neel temper- 
ature of 325 K where the Cu +2 3d 9 ion carries a local 
magnetic moment with spin S = 1/2 due to the sin- 
gle hole in the Cu 3d-shell. 14-17 Then hole doping such 
as in I^-zSr^CuC^ leads to high-T c superconductiv- 
ity where the localized Cu magnetic moment is retained 
and the antiferromagnetic correlations between them are 
likely involved with the superconducting mechanism. 14 
Considering this feature, we speculated that a 122-type 
compound containing Cu +2 local moments with S — 
1/2 completely replacing the Fe in the Fe-sites of the 
122-type compounds might bridge the gap between the 
two families of high-T c superconductors and enrich our 
knowledge of the physics of high-T c superconductivity 
in general. This would complement our recent discov- 
ery that Bai_ a; K a ;Mn2As2 is such an antiferromagnetic 
local-moment metal where the Mn +2 ions have a large 
spinS = 5/2. 27 < 28 

Electronic structure calculations for SrCri2As2 and 
BaCu2As2 by D. J. Singh predicted that the Cu 3d bands 
are narrow and are located about 3 eV below the Fermi 
energy E-p. 29 These calculations also found that there is 
little contribution to the density of states at Ep from 
the Cu 3d orbitals. Therefore these compounds are pre- 
dicted to be sp-band metals with the Cu atoms having a 
formal oxidation state of Cu +1 and a nonmagnetic and 
chemically inert 3d 10 electronic configuration. 29 

To investigate these issues, we synthesized single crys- 
tals of the previously known 30-33 Cu-based compounds 
SrCu2As2, SrCu2Sb2 and BaCu2Sb2 and measured their 
structural, magnetic, thermal and electronic transport 
properties, and report here our results. We also syn- 
thesized and studied crystals of SrCu2(Aso.84Sbo.i6)2- 
To our knowledge, detailed investigations of the phys- 
ical properties of these compounds have not been pre- 
viously carried out. Pfisterer and Nagorsen reported in 
1983 that the \ of a polycrystalline sample of SrCu2As2 
was diamagnctic and nearly independent of T from 80 to 
400 K. 33 

Our measurements confirm the above theoretical 
prediction 29 that SrCu2As2 is an sp-band metal and we 
find that it also applies to the other three compounds 
studied here as well. Thus the Cu ions have a formal 
oxidation state of Cu +1 with a nonmagnetic 3d 10 elec- 
tronic configuration. Therefore these compounds cannot 
be considered to be a bridge to the high-T c cuprates. 




FIG. 1: (Color online) Scanning electron micrograph of a clus- 
ter of intergrown SrCu2(Aso.84Sbo.i6)2 crystals. 



A matter of great current interest is the nature, extent 
and even existence of current carrier doping on substitut- 
ing various transition metals for Fe in the high-T c parent 
compounds (Ca, Sr, Ba)Fe2As 2 because of their relevance 
to understanding the changes in the magnetic and super- 
conducting transition temperatures with concentration of 
the substituting elements. We will give a brief overview of 
this topic and then discuss how our results on SrCu2As2 
relate to the previous results. We suggest that the elec- 
tronic character of the Cu and the strength of the As- As 
interlayer bonds both strongly change between weakly 
Cu-substituted BaFe2As2 and pure BaCu2As 2 , perhaps 
via a first-order lattice instability such as a miscibility 
gap in the Ba(Fe2_ 2; Cu. E )2As 2 system. 

The remainder of this paper is organized as fol- 
lows. The experimental details are described in 
Sec. II. Section III describes the crystallographic stud- 
ies. The physical properties of SrCu2As2, SrCu2Sb 2 , 
SrCu2(Aso.84Sbo.i6)2 and BaCu2Sb 2 are described in 
Sees. IV, V, VI and VII, respectively. A discussion of 
our results in the context of previous work is given in 
Sec. VIII. In Sec. VIII A we discuss the occurrence or not 
of X-X interlayer covalent bonding in (Ca,Sr,Ba)M 2 A2 
compounds with the ThCr 2 Si 2 structure, where M = Cr, 
Mn, Fe, Co, Ni and Cu and X = As or P, and con- 
clude that SrCu 2 As 2 is in the collapsed tetragonal phase 
with short As- As interlayer bonds. In Sec. VIII B we dis- 
cuss the relationships between the X-X bonding and the 
magnetic properties. A brief review of the doping proper- 
ties that occur upon substituting other transition metals 
for Fe in BaFe2As2 is given in Sec. VIII C, together with 
a discussion of how our results on SrCu 2 As 2 impact this 
discussion. A summary and our conclusions are given in 
Sec. IX. 
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TABLE I: The average chemical compositions (atomic %) 
obtained from EDX analyses of SrCu2As2, SrCu2Sb2, 
SrCu2(Aso.84Sbo.i6)2 and BaCu2Sb2 crystals. The system- 
atic error in the EDX measurements is estimated to be ~ 2-3 
atomic %. 



Compound 


Sr or Ba 


Cu 


As 


Sb 


SrCu 2 As 2 


18.3 


41.5 


40.2 




SrCu 2 Sb 2 


19.4 


38.9 




41.8 


SrCu2(Aso.84Sb .16)2 


18.8 


40.9 


34.0 


6.3 


BaCu 2 Sb 2 


19.5 


39.4 




41.1 



II. EXPERIMENTAL DETAILS 

Single crystals of SrCu 2 As 2 , SrCu 2 Sb 2 , 
SrCu 2 (As .84Sbo.i6)2 and BaCu 2 Sb 2 were grown by 
the self-flux growth technique using the high purity 
elements Sr (99.95%) and Ba (99.99%) from Sigma 
Aldrich, and Cu (99.999%), As (99.99999%) and Sb 
(99.999%) from Alfa Acsar. Prcreacted CuAs or CuSb 
were used as flux. For the growth of SrCu 2 As 2 and 
SrCu 2 Sb 2 crystals, the Sr and flux (CuAs or CuSb) 
were taken in a 1:5 molar ratio, placed in alumina 
crucibles and sealed inside evacuated quartz tubes. The 
crystal growth was carried out by heating the samples 
to 1100 °C at a rate of 60 °C/h, holding there for 
12 h and then cooling to 800 °C at a rate of 2.5 °C/h. 
The crystals were separated from the flux by decanting 
the flux with a centrifuge at 800 °C, yielding shiny 
plate-like crystals. While the sizes of the SrCu 2 Sb 2 
crystals were typically 3 x 2.5 x 0.3 mm 3 , one of the 
SrCu 2 As 2 crystals was larger (6x3x1 mm 3 ). In 
our attempt to grow crystals of the mixed compound 
SrCu 2 AsSb starting with Sr:CuAs:CuSb in a 1:2.5:2.5 
molar ratio with the same heating profile as above, 
wc instead obtained SrCu 2 (Aso.84Sbo.i6)2 crystals as 
determined from energy-dispersive x-ray (EDX) analysis. 
In this case scanning electron microscope (SEM) images 
revealed intergrowths of shiny plate-like crystals stacked 
on top of each other (Fig. 1). The physical properties of 
this composition were measured on such an assembly of 
stacked crystals. 

For the growth of BaCu 2 Sb 2 crystals, Ba and CuSb 
were taken in a 1:5 molar ratio and heated to 1000 °C at 
a rate of 80 °C/h, held there for 15 h, cooled to 700 °C 
at a rate of 2.5 °C/h, and then the flux was decanted at 
700 °C to obtain crystals of typical size 2.5 x 2 x 0.5 mm 3 . 
We did not succeed in synthesizing single crystals of 
BaCu 2 As 2 using the self-flux growth or with Sn as a flux 
with similar heating profiles. 

The chemical compositions of the crystals were deter- 
mined using a JEOL SEM equipped with an EDX ana- 
lyzer, which showed that these are close to the respec- 
tive stoichiometric compositions SrCu 2 As 2 , SrCu 2 Sb 2 , 
SrCu 2 (Aso.84Sbo.i6)2 and BaCu 2 Sb 2 . The average chem- 
ical compositions obtained from the EDX analyses of 
the single crystals are presented in Table I. The crys- 



tal structures of the samples were determined by powder 
x-ray diffraction (XRD) measurements that were carried 
out on powdered crystals using a Rigaku Gcigcrflcx x-ray 
diffractometer and Cu K Q radiation. 

The specific heat measurements were performed by the 
relaxation method in zero applied magnetic field in a 
Quantum Design physical properties measurement sys- 
tem (PPMS). The electrical resistivity was measured in 
the a6-plane of each crystal in zero applied magnetic field 
by the standard four-probe ac technique using the ac 
transport option of the PPMS. The 50 /xm diameter plat- 
inum wire electrical leads were attached to the crystals 
using silver epoxy. 

Magnetization and magnetic susceptibility measure- 
ments were performed using a Quantum Design super- 
conducting quantum interference device magnetic prop- 
erties measurement system (MPMS). The single crystals 
were mounted either on a small piece of plastic or on a 
thin quartz rod using a small amount of GE 7031 var- 
nish. We separately measured the magnetic moment of 
the sample holder (quartz rod/plastic and GE varnish) 
and subtracted it from the measured magnetic moment 
to obtain the sample contribution. 



III. CRYSTALLOGRAPHY 

The XRD data on crushed single crystals of SrCu 2 As 2 , 
SrCu 2 Sb 2 , SrCu 2 (Aso.84Sbo.i6)2 and BaCu 2 Sb 2 were an- 
alyzed by Rietveld refinement using the FullProf 34 
software. The room-temperature XRD patterns along 
with the Rietveld fit profiles are presented in Fig. 2 for 
SrCu 2 As 2 , SrCu 2 Sb 2 and SrCu 2 (Aso.84Sbo.i6)2 and in 
Fig. 3 for BaCu 2 Sb 2 . The XRD data indicate that the 
crushed crystals are single-phase. The unindexed peaks 
marked with stars arise from small amounts of flux at- 
tached to the surfaces of the crystals. The Rietveld re- 
finements confirmed that SrCu 2 As 2 crystallizes in the 
ThCr 2 Si 2 -type body-centered tetragonal (bet) structure 
(space group IA/mmm) and SrCu 2 Sb 2 crystallizes in 
the CaBe 2 Ge 2 -type primitive tetragonal structure (space 
group PA/nmm). As seen in Fig. 4(a), the structure of 
SrCu 2 As 2 contains stacked square lattices of Cu atoms, 
as in the high T c cuprates, but the coordination of cop- 
per by the anions is quite different in the two types of 
material. On the other hand, the structure of SrCu 2 Sb 2 
contains two distinct types of Cu square lattice with dif- 
ferent lattice parameters that are rotated by 45° from 
each other and stacked in a specific sequence, as seen in 
Fig. 4(b). 

Both the ThCr 2 Si 2 -type and CaBe 2 Ge 2 -type struc- 
tures are ternary derivatives of the BaAl 4 structure 35 
and consist of layers of A (Th, Ca), T (Cr, Be) and 
X (Si, Ge) atoms stacked along the crystallographic c- 
axis. While the A atoms form a bet sublattice in both 
cases, the two structures differ in the arrangement of lay- 
ers of T and X atoms as can be seen from Figs. 4(a) 
and 4(b). In the ThCr 2 Si 2 -type structure the layers of 
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TABLE II: Crystallographic and Rietveld refinement parameters obtained from powder XRD data for crushed crystals of 
SrCu2As2, SrCu 2 Sb 2 , SrCu2(Aso.84Sb .i6)2 and BaCu 2 Sb 2 . If no value for the lattice parameter b is listed, then b — a. Error 
bars for the last digit of a quantity are given in parentheses and literature references are given in square brackets. 





SrCu 2 As2 


SrCu 2 Sb 2 


SrCu2 ( Aso. 84 Sbo. i6 ) 2 


BaCu 2 Sb 2 


BaCu 2 Sb 


Structure 


ThCr 2 Si 2 -type 


CaBe2Ge2-type 


ThCr 2 Si 2 -type 








tetragonal 


tetragonal 


tetragonal 


tetragonal 


orthorhon 


Space group 


IA/mmm 


PA/nmm 


IA/mmm 


IA/mmm 


Immm 


Lattice parameters 












a (A) 


4.2725(1) 
4.279(1) [31] 
4.271(1) [30] 


4.5162(2) 
4.51(1) [32] 


4.2940(1) 


4.6438(5) 
4.655(1) [31] 


4.6224(5) 


6(A) 










4.6679(5) 


C (A) 


10.2000(3) 
10.215(1) [31] 


10.9008(5) 
10.91(2) [32] 


10.3190(4) 


32.600(4) 
32.709(6) [31] 


32.606(4) 


Vcell (A 3 ) 


10.18(2) [30] 










186.19(1) 


222.33(2) 


190.26(1) 


703.0(1) 


703.5(1) 


Refinement quality 












x 2 


9.72 


9.13 


12.3 


4.98 


4.00 


R P (%) 


3.70 


6.81 


4.59 


5.13 


4.77 


^?wp (%) 


5.71 


9.52 


7.32 


6.98 


6.26 



TABLE III: Atomic coordinates obtained from the Rietveld 
refinements of powder XRD data for powdered crystals of 
SrCu2As2, SrCu2Sb2, SrCu2(Aso.84Sbo.i6)2 and BaCu2Sb2. 



Atom 



Wyckoff x y 



symbol 


SrCu2As2 (IA/mmm) 










Sr 


2a 











Cu 


4d 





1/2 


1/4 


As 


4e 








0.3789(1) 


SrCu 2 Sb 2 (PA/nmm) 










Sr 


2c 


1/4 


1/4 


0.2353(5) 


Cul 


2a 


3/4 


1/4 





Cu2 


2c 


1/4 


1/4 


0.6330(8) 


Sbl 


2b 


3/4 


1/4 


1/2 


Sb2 


2c 


1/4 


1/4 


0.8707(4) 


SrCU2(Aso.84Sb .16)2 










Sr 


2a 











Cu 


4d 





1/2 


1/4 


As/Sb 


4e 








0.3781(2) 


BaCu 2 Sb 2 (IA/mmm) 










Bal 


2a 











Ba2 


4e 








0.6654(4) 


Cul 


4e 








0.7947(7) 


Cu2 


8g 





1/2 


0.0830(6) 


Sbl 


4d 





1/2 


1/4 


Sb2 


4e 








0.1234(3) 


Sb3 


4e 








0.5428(3) 


BaCu2Sb2 (Immm) 










Bal 


2a 











Ba2 


4i 








0.6650(3) 


Cul 


4i 








0.7933(6) 


Cu2 


4j 


1/2 





0.0808(9) 


Cu3 


4j 


1/2 





0.5880(8) 


Sbl 


4j 


1/2 





0.7484(8) 


Sb2 


4i 








0.1230(3) 


Sb3 


4i 








0.5426(3) 



atoms are in the order Th-Si-Cr-Si-Th-Si-Cr-Si-Th along 
the c-axis, whereas in the CaBe2Ge2-type structure the 
order of atomic layers is Ca-Ge-Be-Ge-Ca-Be-Ge-Be-Ca. 
Thus in between the consecutive Ca-planes the Ge-Be-Ge 
slab changes into a Be-Ge-Be slab which leads to a loss 
of mirror symmetry about the A-plane in the CaBe2Ge2- 
type structure that is present in the ThCr 2 Si2-type struc- 
ture. SrCu2(Aso.84Sbo.i6)2 is found to crystallize in the 
bet ThCr2Si2-type structure. 

On the other hand the XRD data for BaCu2Sb2 
could not be described by either a ThCr2Si2-type or 
a CaBe2Ge2-type structure. Instead, as reported by 
Diinner et al., 31 the structure of BaCu2Sb2 is an or- 
dered intergrowth of ThCr 2 Si2- and CaBe2Ge2-types of 
unit cells. The unit cell of BaCu2Sb2 consists of three 
of these unit cell blocks stacked along the c-axis, as 
shown in Fig. 4(c), leading to a large c = 32.6 A, 
with a ThCr2Si2-type block sandwiched between two 
CaBe2Ge2-type blocks resulting in an overall symme- 
try of IA/mmm. A Rietveld refinement of the XRD 
data of BaCu2Sb2 with this structural model is shown in 
Fig. 3(a). The overall agreement between the experimen- 
tal XRD data and the calculated pattern is reasonable. 
However, we observe splittings of some of the XRD peaks, 
as illustrated in the inset of Fig. 3(a), which suggest that 
the symmetry is lower than the reported one. 

We attempted to refine the XRD data for BaCu 2 Sb 2 
with a lower orthorhombic (Immm) symmetry, which 
is a subgroup of the IA/mmm space group. The Ri- 
etveld refinement profile with this structural model is 
shown in Fig. 3(b). In terms of the goodness of fit, this 
structural model proved better as deduced from the re- 
spective x 2 values in Table II, but it could not account 
for the magnitudes of the splittings of the peaks, as il- 
lustrated in the inset of Fig. 3(b). Further lowering of 
the symmetry to monoclinic space group C2/m, which 
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FIG. 2: (Color online) Powder x-ray diffraction patterns of 
(a) SrCu 2 As 2 , (b) SrCu 2 Sb 2 and (c) SrCu 2 (As .84Sb .i6)2 
recorded at room temperature. In (a) and (c), the solid 
lines through the experimental points are the Rietveld refine- 
ment profiles calculated for the ThCr 2 Si 2 -type body-centered 
tetragonal structure (space group 74/mmm), and in (b) the 
solid line through the data points is the Rietveld refinement 
profile for the CaBe 2 Ge 2 -type primitive tetragonal structure 
(space group P4/nmm). In (a), (b) and (c), the short ver- 
tical bars mark the fitted Bragg peak positions. The lower- 
most curves represent the differences between the experimen- 
tal and calculated intensities. The unindexed peaks marked 
with stars correspond to peaks from the flux. 
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FIG. 3: (Color online) Powder x-ray diffraction pattern of 
BaCu 2 Sb 2 recorded at room temperature, (a) The solid line 
through the experimental points is the Rietveld refinement 
profile calculated for the tetragonal structure (space group 
I4/mmm). (b) Rietveld refinement profile calculated for an 
orthorhombic structure (space group Immm). (c) Le Bail 
profile fit for a monoclinic structure (space group P2\/c). In 
(a), (b) and (c), the short vertical bars mark the fitted Bragg 
peak positions, the lowermost curves represent the differences 
between the experimental and calculated intensities, and the 
insets show a small section of the XRD pattern highlighting 
an exemplary splitting and fitting of one of the peaks. 
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FIG. 4: (Color online) (a) ThCr2Si2-type body-centered 
tetragonal crystal structure (74/mmm) of SrCu2As2. 
(b) CaBe2Ge2-type primitive tetragonal crystal structure 
(PA/nmm) of SrCu2Sb2- The order of the Cu and Sb layers 
in the lower half of the unit cell is reversed in (b) compared 
with (a). In order to compare the structures, the origin of 
the SrCu 2 Sb 2 unit cell is shifted by (1/4, 1/4, 1/4) from that 
for the P4:/nmm space group, (c) Body-centered tetragonal 
representation of BaCu2Sb2 (space group IA/mmm), consist- 
ing of a central ThCr2Si2-type block sandwiched between two 
CaBe2Ge2-type blocks, resulting in a large c-axis parameter 
c = 32.586 A (see Table II). 



is a subgroup of Immm, gave a better refinement qual- 
ity but could not account for all of the peaks present. 
Further lowering of the symmetry to space group P2i/c, 
which is a subgroup of C2/m, allowed us to fit all the 
peak positions as illustrated in the inset of Fig. 3(c). A 
Le Bail profile fit for this P2\/c monoclinic structure 
with a = 4.6414(3), b = 4.6418(3), c = 32.680(2) A and 
= 90.616° is shown in Fig. 3(c). A similar lowering 
of symmetry to P2\/c was necessary for the 122-type 
compound SrRli2As2 which exhibits polymorphism from 




T(K) 

FIG. 5: (Color online) Heat capacity C p of a SrCu2As2 sin- 
gle crystal versus temperature T. The solid red curve is a 
fit of the data by the sum of the contributions from the De- 
bye lattice heat capacity CvDcbye(T) and predetermined elec- 
tronic heat capacity 7T according to Eq. (12). Inset: C p /T 
vs. T 2 below 5 K. The straight line is a fit of the data for 
1.85 K < T < 3.5 K by Eq. (1). 



the high-T tetragonal ThCr 2 Si2-type structure to a low-T 
P1\jc monoclinic structure below 350 °C. 36 We conclude 
that the most likely space group is the monoclinic P2±/c 
that we have identified. The detailed crystal structure 
(determination of the precise atomic positions) remains 
to be determined. Our attempt to solve the structure 
using single crystal x-ray diffraction failed because the 
crystal quality was insufficient to obtain reliable data. 



The crystallographic and refinement parameters of all 
four compounds are listed in Tables II and III. The lat- 
tice parameters obtained from our refinements are in rea- 
sonable agreement with available literature values, 30-32 
which are also listed in Table II. During the refinements 
the thermal parameters B were kept fixed to B = 
since the lattice parameters and z parameters were in- 
sensitive to change in B within the error bars. Further- 
more, the occupancies of atomic positions were kept fixed 
at the value of unity as there were no improvements in 
the qualities of fit on varying the occupancies by small 
amounts from unity. The z coordinate of the As atoms in 
SrCu2As2, zas — 0.3789(1), is very close to the reported 
values z As = 0.3785(1) (Ref. 31) and z As = 0.377. 33 
The ^-coordinates z$r = 0.2353(5), zq U 2 = 0.6330(8) and 
zsb2 = 0.8708(4) of SrCu2Sb2 are also very close to the 
reported values zs r = 0.2369(6), zq U 2 = 0.6380(6) and 
z S b2 = 0.8695(6). 32 



7 



IV. PHYSICAL PROPERTIES OF SrCu 2 As 2 
CRYSTALS 

A. Heat Capacity 

The heat capacity C p of SrCu2As2 as a function of T 
measured at constant pressure p is presented in Fig. 5. 
No evidence for any phase transitions was observed down 
to 1.85 K. At 300 K the heat capacity attains a value 
of 123.5 J/molK which is close to the expected classical 
high-T Dulong-Petit lattice heat capacity value Cy = 
3nR = 15R = 124.7 J/molK at constant volume V, 37 ' 38 
where R is the molar gas constant and n = 5 is the 
number of atoms per formula unit (f.u.) of SrCu2As 2 . 
The inset in Fig. 5 shows the low-T data plotted as C p /T 
versus T 2 allowing a conventional fit by 



C P (T) 



(1) 



where 7 is the Sommcrfeld electronic linear specific heat 
coefficient and f3 is the coefficient of the Debye T 3 lattice 
heat capacity at low T. A linear fit of the data below 
3.5 K according to Eq. (1), shown as the straight line in 
the inset of Fig. 5, gives the parameters 



7 = 2.2(2) mJ/molK 2 
= 0.47(3) mJ/molK 4 . 



(2) 



We obtained the value of the Debye temperature 0d from 
/3 using the relation 37 



127r 4 A^fc B n \ 1/3 



(3) 



where Na is Avogadro's number and fcs is Boltzmann's 
constant, yielding 



9 D = 274(6) K. 



(4) 



The Sommcrfeld coefficient 7 can be used to estimate 
the density of states at the Fermi level V{E F ) for both 
spin directions according to 37 
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kD(£ F )(l + A e -ph), 



(5) 



where A e - P h is the electron-phonon coupling constant 
which we set to zero and obtain a value 



V(E F ) = 0.93 states/eVf.u. 



(6) 



for both spin directions. This value of T>(E F ) is less 
than the value 1.53 states/eVf.u. for both spin direc- 
tions predicted from the band structure calculations for 
SrCu 2 As 2 . 29 

The V(E F ) can be used to estimate the band effective 
mass to* using the relation 



V(E F ) = iV band 



*Vf.u.k F 



(7) 



where T>(E F ) is in units of states/erg f.u. with the quan- 
tities on the right-hand side in cgs units, k F is the Fermi 
wave vector, Vf. u . is the volume per formula unit and 
Aband is the number of degenerate bands at E F with the 
same to* and Fermi velocity v F . The quantity in paren- 
theses on the right-hand side corresponds to the single- 
band relation. 37 Further, using k F — m*v F /h, Eq. (7) 
can be rearranged to yield 



h\ 2 



V(E F ) 



Vf.u. TO 2 V F Aband 



(8) 



where m e is the free electron mass. Using the average 
value 



(v F ) = 2.7 x 10 7 cm/s 



(9) 



from the band structure calculation 29 and Vf. u . = V ce \\/Z 
where Z = 2 is the number of formula units per unit cell, 
we obtain 



TO 
TO„ 



= 1.86 



V(E F 



1 1/2 



band 



(10) 



where V{E F ) is now in units of states/eVf.u. for both 
spin directions. The band structure calculations show 
that Aband = 2 for SrCu2As2, but there are three distinct 
Fermi surfaces arising from the two bands. 29 Hence tak- 
ing Aband = 2 or 3 and using V(E F ) = 0.93 states/eVf.u. 
for both spin directions from Eq. (6), the effective mass 
estimated from Eq. (10) is 



TO* = 1.27 TO C (Aband = 2) 

TO* = 1.04to c (Aband = 3). 



(11) 



The ratio to*/to c ~ 1 is consistent with the prediction 29 
that SrCu 2 As 2 is an sp-band metal. 

The heat capacity data of a nonmagnetic metal over 
an extended T range can be approximated by 



C P (T) = 7T + nCV Dcbyc(7 1 ), 



(12) 



where 7 is fixed to be the Sommerfeld electronic specific 
heat coefficient found above, Cy Debye CO is the Debye 
lattice heat capacity due to acoustic phonons at con- 
stant volume V per mole of atoms, and n is the num- 
ber of atoms per f.u. (n = 5 here). The Debye function 
CVDcbyc(T) is given by 38 



CVDcbyc(T) = 9i? ( — 



T \ 3 f eD/T x 4 e x 



(e* - l) 2 



dx. (13) 



The heat capacity data in Fig. 5 were least-squares 
fitted by Eqs. (12) and (13) over the T range 2-300 K 
using the Pade approximant function for Cy Debye (T) de- 
veloped by us in Ref. 39, with Gd being the only ad- 
justable parameter. The solid curve in Fig. 5 represents 
the best fit of the experimental data in Fig. 5 by Eqs. (12) 
and (13), yielding 6 D = 246(1) K. This value of 9 D is 
different from the value of 8d = 274(6) K obtained in 
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TABLE IV: The linear specific heat coefficients 7 and the coefficient j3 of the T 3 term in the low-T heat capacity, and the density 
of states at the Fermi energy T>(Ef) for both spin directions for SrCu2As2, SrCu2Sb2, SrCu2(Aso.84Sbo.i6)2 and BaCu2Sb2 
single crystals. The Debye temperatures <do obtained at low T and for all T from heat capacity measurements and the Debye 
temperature 0r obtained from fitting electrical resistivity data, respectively, are also listed. 



Compound 


7 


P 


V(E F ) 


e D (K) 


e D (K) 


e R (K) 




(mJ/molK 2 ) 


(mJ/molK 4 ) 


(states/eV f.u.) 


from low-T 


from all T 


from p(T) 


SrCu 2 As 2 


2.2(2) 


0.47(3) 


0.93 


274(6) 


246(1) 


245(4) 


SrCu 2 Sb 2 


3.9(2) 


0.76(3) 


1.65 


234(3) 


214(1) 


153(1) 


SrCu 2 (Aso.84Sb .i6)2 


2.3(2) 


0.83(2) 


0.97 


227(2) 


237(1) 


225(6) 


BaCu 2 Sb 2 


3.5(2) 


0.65(2) 


1.49 


246(3) 


204(1) 


160(2) 



Eq. (4) from the coefficient (3 deduced from the low-T 
heat capacity data. The difference between the two val- 
ues may arise due to the approximations made in Debye 
theory that result in a T-dependent 9d, so the fit by 
Eq. (12) gives a value of 6d that is averaged over T. A 
contribution to the difference also comes from the fact 
that the heat capacity was measured at constant pres- 
sure whereas the Debye theory describes the lattice heat 
capacity at constant volume. The parameters obtained 
from our heat capacity fits are summarized in Table IV. 
Also listed is the Debye temperature 0r obtained from 
fitting our experimental p(T) data for SrCu2As2 by the 
Bloch-Gruneisen model in Sec. IV C below. 



B. Magnetization and Magnetic Susceptibility 

The zero-field-cooled (ZFC) magnetic susceptibilities 
\ = M /H of a SrCu2As2 single crystal as a function 
of T from 1.8 to 340 K in an applied magnetic field 
H = 3.0 T aligned along the c-axis (\c H || c) and in 
the a6-plane (xab, H _L c) are presented in Fig. 6, where 
M is the measured magnetization of the crystal. The x 
has a negative sign and is anisotropic with Xab > Xc over 
the whole T range. This kind of anisotropy is common in 
the iron arsenide family; e.g., both BaFe2As2 (Ref. 40) 
and SrFe2As2 (Ref. 41) exhibit Xab > Xc- The pow- 
der and temperature average of our anisotropic x va l~ 
ues in Fig. 6 is significantly more diamagnetic (negative) 
than the T-independent value between 80 and 400 K of 
X = —5.4 x 10~ 5 cm 3 /mol reported previously for a poly- 
crystalline sample. 33 This difference may be due to a sat- 
uration magnetization of ferromagnetic (FM) impurities 
in the latter sample that was incompletely corrected for. 

Upturns are observed in the x(T) data in Fig. 6 at low 
T for both Xab and % c that are likely due to the pres- 
ence of a small amount of saturable paramagnetic (PM) 
impurities as deduced from the M(H) isotherm data dis- 
cussed below. Since SrCu2As2 is weakly diamagnetic, 
even trace amounts of magnetic impurities can make sig- 
nificant contributions to the x an d M data. Therefore, a 
correction for the magnetic impurity contribution is nec- 
essary to determine the intrinsic magnetic behavior of 
the compound. 

The isothermal M data for SrCu2As2 as a function of 
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FIG. 6: (Color online) Zero-field-cooled magnetic susceptibil- 
ity x of a SrCu2As2 single crystal versus temperature T in 
the temperature range 1.8 - 350 K measured in a magnetic 
field of 3.0 T applied along the c-axis (xc, H || c) and in 
the afc-plane (Xab, H 1 c) (solid symbols). The open sym- 
bols represent the intrinsic susceptibility of SrCu2As2 after 
correcting for the ferromagnetic and paramagnetic impurity 
contributions as described in the text. 



H measured at eight temperatures between 1.8 and 300 K 
with H applied along the c-axis (M c , H || c) and in the 
a6-plane (M a t, Hie) are presented in Fig. 7. The 
M(H) data are diamagnetic and reveal anisotropic be- 
havior with M a b(H) > M C (H), which is consistent with 
the above observation Xab > Xc- The presence of slight 
nonlinearities in the M(H) isotherms at low fields, par- 
ticularly at low T, reveals the presence of saturable fer- 
romagnetic FM and/or paramagnetic PM impurities in 
the sample. 

Since the M(H) contribution of a FM impurity is typ- 
ically nonlinear only at low fields H < 1-2 T one can 
estimate the FM impurity saturation magnetization M s 
from the high-field M (H) data. Furthermore, typical FM 
impurities have Curie temperatures significantly above 
room temperature, and saturable PM impurities do not 
show nonlinear M(H) behaviors above ~ 25 K. We 
therefore determined the M s of the FM impurities by 
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FIG. 7: (Color online) Isothermal magnetization M of a 
SrCu2As2 single crystal versus magnetic field H measured at 
the indicated temperatures for H applied (a) perpendicular to 
the c-axis (M a b, H 1 c) and (b) along the c-axis (M c , H || c). 
The solid curves are fits by Eq. (15). 



fitting the M(H) data at H > 2.0 T for T > 25 K by 

M(H) = M s + X H, (14) 

where \ is the susceptibility of the sample for a particu- 
lar field orientation that can still contain the contribution 
from paramagnetic PM impurities as discussed next. The 
M s values obtained from the fits were found to be inde- 
pendent of T for T > 25 K to within the error bars. The 
values are listed in Table V. The M s c value corresponds 
to the saturation magnetization of only 4 molar ppm of 
Fe metal impurities. 

The PM impurity contribution to M can be estimated 
by fitting the M(H) data for each field direction by 

M(T, H) = M s + X oH + / imp A/ SimpJ B Sim >), (15) 

where M s is the FM impurity saturation magnetization 
value found above, xo is the intrinsic susceptibility of 



TABLE V: Parameters obtained for SrCu2As2 and SrCu2Sb2 
from fitting M(H) isotherms at 1.8 K by Eqs. (14) and (15), 



where 9i„ 



0. Here M s is the saturation magnetization 



of ferromagnetic impurities, \o is the intrinsic susceptibility, 
and /i m p and Si mp are the molar fraction and spin of the 
paramagnetic impurities, respectively. 



Compound 


field 


M s 


Xo 


fi mp S j mp 




direction 


( G cm 3 \ 
I mnl ) 


(10- 5 sei) 


(io- 5 ) 


SrCu 2 As 2 


H _L c 


0.00(2) 


-7.62(2) 


3.05(8) 




H\\c 


0.05(3) 


-9.32(4) 


4.1(2) 


SrCu 2 Sb 2 


H _L c 


0.00(5) 


-17.1(1) 


10.9(3) 




H || c 


0.02(3) 


-5.56(3) 


2.2(1) 



the compound, /i mp is the molar fraction of PM impuri- 
ties, M Si = A^AffimpMB'S'imp is the PM impurity satura- 
tion magnetization, is Avogadro's number, /ie is the 
Bohr magneton, and gi mp and £i mp are the spectroscopic 
splitting factor (g-factor) and the spin of the impurities, 
respectively. The Brillouin function Bs ilnp for the PM 
impurities is given by 



B 8bap (x) 



— — \ (2S imp + 1) coth \(2S imp + 1)| 



where 



coth 



_ gimp^BH 




(16) 



(17) 



k B (T - imp ) 

A Weiss temperature 9 lmp is included in the Brillouin 
function in order to take weak interactions between PM 
impurity spins into account in a mean-field way. In par- 
ticular, when Bs ilap (x) is expanded to first order in x for 
x <C 1, one obtains a Curie- Weiss law for the low-field 
susceptibility of the impurity spins with Weiss tempera- 
ture #i mp . In order to reduce the number of fitting pa- 
rameters, the impurity (/-factor was set to g- lmp = 2 for 
this compound as well as for SrCu2Sb2 in Sec. VB below. 

The solid curves in Fig. 7 are the fits of the M(H) data 
by Eq. (15) that were obtained using the above M s ° & and 
Af s c values for H _L c and H || c, respectively. Since M s ah 
is zero we fitted the M(H) data for H _L c by Eq. (15) 
in the whole range < H < 5.5 T. However, since M s c is 
nonzero we fitted the M(H) data for H || c only in the 
range 1.0 < H < 5.5 T. We found that setting 6>j mp = 
gives a good fit for both field directions, so in the final 
fits we set flimp = 0. 

The parameter xo an d the product of the parame- 
ters /i m p and iSi mp obtained from the fits of the M(H) 
isotherms at 1.8 K for H || c and H _L c in Fig. 7 by 
Eq. (15) are listed in Table V. From Table V, both M s 
and the product /impShnp are different for H || c com- 
pared with Hie, which indicates that the magnetiza- 
tions of the FM and PM impurities are both anisotropic. 
The origins of these anisotropies are not clear. 
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From Eq. (15), the magnetic impurity contribution to 
the magnetization, for fields above the saturation field of 
about 1 T for the FM impurities, is 

M imp (iJ,T) = M s + / imp M Simp B Simp (x). 

The x(T) = M(T)/H data in Fig. 6 were taken with 
H = 3 T. Therefore, this field was above the saturation 
field of the FM impurities and one can correct for the 
contributions of both the FM and PM impurities to ob- 
tain the intrinsic susceptibility versus T according to 



xoCO = 



M(T) - M lmp (H,T) 
H 



In this way we obtained the intrinsic Xo{T) data shown 
as the open symbols in Fig. 6. The magnetic impurities 
are seen to have little influence on the measured suscep- 
tibility except below about 50 K where the PM impurity 
contribution to M(T) becomes significant. 

The largest value of /impSimp = 1.09 x 10~ 4 in Table V, 
which is for Xab of SrCu 2 Sb 2 (see below), corresponds to 
a very small concentration of paramagnetic impurities. 
For the worst case assumption <Si m p = 1/2, the impu- 
rity concentration is /i mp = 2.2 x 10~ 4 = 0.022 mol% = 
220 molar ppm. If Si mp = 5/2, one obtains the five-times 
smaller value / imp = 4.4 x 10~ 5 = 0.0044 mol% = 44 mo- 
lar ppm. The values of M s < 0.05 Gcm 3 /mol in Table V 
correspond to M s < 9.0 x 10~ 6 // B /f.u. Thus if S = 1/2 
and g = 2, this corresponds to a ferromagnetic molar im- 
purity concentration of < 9.0 molar ppm. On the other 
hand, if S = 5/2 and g = 2, this corresponds to a fer- 
romagnetic molar impurity concentration of < 1.8 molar 
ppm. As stated in Sec. II, the purity of the Sr used in 
the synthesis of the SrCu 2 As 2 and SrCu 2 Sb2 crystals was 
99.95 mol% (metals basis), i.e., containing 0.05 mol% of 
elemental metal impurities. Thus our values of /impSimp 
and M s in Table V are consistent with (less than) the ele- 
mental metal impurity level supplied by the source of our 
Sr starting material. The nature of these magnetic impu- 
rities is unknown, and this information was not provided 
by the supplier of the Sr starting material. 

The residual low-T upturns in Xab and Xc in Fig. 6 
(and in Fig. 10 below) after correction for the contribu- 
tion of the paramagnetic impurities are very small and 
are most likely due to imperfect subtractions of the para- 
magnetic impurity contributions. The small changes in 
the magnitudes of x between 300 and 350 K in Fig. 6 (and 
in Figs. 10 and 14 below) are believed to be experimen- 
tal artifacts of our SQUID magnetometer. Such small 
changes are seen in many samples of different types that 
we have measured on this instrument in this T range. 
The origin of this systematic error is unknown. 

The contributions to xo of a nonmagnetic metal are 



XO = Xcorc + XVV + XL + XP, 



(18) 



where the first three terms comprise the orbital suscep- 
tibility contributions. Here Xcore is the diamagnetic sus- 
ceptibility due to the localized core electrons, xvv is the 



TABLE VI: Estimated contributions to the intrinsic angle- 
and temperature-averaged magnetic susceptibilities (xo) of 
SrCu2As2, SrCu2Sb2, SrCu2(Aso.s4Sbo.i6)2 and BaCu2Sb2 
crystals in units of 10 cm /mol. Here xp is the Pauli spin 
susceptibility of the conduction carriers, and the orbital sus- 
ceptibility contributions are the diamagnetism Xcorc of the 
atomic electron cores, the Landau diamagnetism xl of the 
conduction carriers and the Van Vleck paramagnetism xvv- 



Compound 


<Xo) 


Xcorc 


Xp 


XL 


Xvv 


SrCu 2 As 2 


-8.25 


-17.50 


3.00 


-0.62 


6.87 


SrCu 2 Sb2 


-15.58 


-20.79 


5.33 


-1.78 


1.66 


SrCu2(Aso.84Sbo.ie)2 


-6.95 


-18.03 


3.14 


-1.05 


8.99 


BaCu 2 Sb 2 


-8.02 


-23.10 


4.81 


-1.60 


11.9 



paramagnetic Van Vleck susceptibility, and xl is the dia- 
magnetic Landau susceptibility of the conduction carri- 
ers. The last term xp is the paramagnetic Pauli spin 
susceptibility of the conduction carriers. 

The powder- and temperature-average (from 1.8 to 
350 K) of xo is (xo) = VXab + Xc)/3 = -8.25 x 
10~ 5 cm 3 /mol. Since our compounds are closer to be- 
ing covalent metals than ionic metals, we used atomic 
diamagnetic values 42 for Xcorc to obtain Xcorc = —1-75 x 
10~ 4 cm 3 /mol. The xp can be estimated from T>(Ep) 
using the relation 43 



Xp = »%V{E ¥ ). 



(19) 



We then obtain xp = 3.00 x 10~ 5 cm 3 /mol using 
T>(Ep) = 0.98 states/eVf.u. for both spin directions from 
Eq. (6). The xl can be estimated from the relation 43,44 



XL 



I ( m e\ 2 

"3 \m*J 



Xp- 



(20) 



Using m* — 1.27 m c from Eq. (11), one obtains xl = 
—0.62 x 10~ 5 cm 3 /mol. One can then estimate the 
angle- and temperature-averaged xvv using Eq. (18) 
and the above values of Xcorc, Xp an d XL, which give 
Xvv = 6.87 x 10 -5 cm 3 /mol, a physically reasonable 
value. The above contributions to the average intrin- 
sic magnetic susceptibility of SrCu 2 As 2 are summarized 
in Table VI. 



C. Electrical Resistivity Measurements and the 
Bloch-Griineisen Model 

The in-plane electrical resistivity p data for SrCu 2 As 2 
as a function of T are presented in Fig. 8. Metallic char- 
acter is indicated from the positive temperature coef- 
ficient of p(T). The relatively low value (for 122-typc 
compounds) of the residual resistivity 

Po = p(T 0) = 9.3 pn cm 



and a reasonably large residual resistivity ratio 

6 
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FIG. 8: (Color online) In-plane electrical resistivity p of a 
SrCu2As2 crystal versus temperature T. The red solid curve 
is the fit by the Bloch-Griineisen model. 



indicate that our crystal is of good quality. 

We fitted the p(T) data of SrCu 2 As 2 by the Bloch- 
Griineisen (BG) model that describes the electrical re- 
sistivity pbg(T) due to scattering of the conduction elec- 
trons by acoustic lattice vibrations in a monatomic metal 
according to 45 



Pbg (t) = 4^(e R ) 



T 

Or" 



Or/T 



(e x - 1)(1 - e~ x ) 



-dx, 



(21) 

where R is the Debye temperature obtained from a fit 
of experimental resistivity data by the BG theory and 
the prefactor function 7\L(9 R ) is given by 



^ 3 (37T 2 ) 1/3 ft 2 

4n^aMfc B e R 



(22) 



Here Ti is Planck's constant divided by 27r, e is the funda- 
mental electric charge, n co n is the number of conduction 
(valence) electrons per unit cell of volume a 3 containing 
one atom, M is the mass of the one atom per unit cell and 
k-Q is Boltzmann's constant. The prefactor in Eq. (22) is 
h/e 2 = 4108.24 O, in SI units. Thus 7£(6 R ) can be found 
in the favored units of cm by expressing the quanti- 
ties in square brackets in cgs units. The expression for 
TZ(Q R ) in Eq. (22) needs to be modified slightly accord- 
ing to Eq. (26) below in order to apply the BG theory to 
calculating the magnitude of the resistivity of polyatomic 
metals. 

Even for simple s- or sp-metals it is usually not possible 
to accurately fit both the magnitude and T dependence 
of p by the BG model with R as the only adjustable 
parameter. 45 To fit the T dependence of p using the BG 
model we first calculate the BG prediction for p(T = <d R ) 
from Eq. (21) as 39 



TABLE VII: Parameters associated with Bloch-Griineisen fits 
to the resistivities p within the ab-plane of the listed single 
crystals. Here po is the residual resistivity for T — Y 0, 0r is 
the Debye temperature determined from resistivity measure- 
ments, p(Qr) is the fitted value of p at T = Qn, and 1Z(Or) 
is the quantitiy defined in Eq. (22) and is obtained from the 
fitted value of p(Qr) using Eq. (23). 



Compound 


Po 


e R 


p(e R ) 


ft(6 R ) 




(pQ cm) 


(K) 


(pSl cm) 


(pQ cm) 


SrCu 2 As 2 


9.90(7) 


245(4) 


34.9(6) 


36.9 


SrCu 2 Sb 2 


13.64(1) 


153(1) 


13.8(1) 


14.6 


SrCU2(Aso.84Sbo.l6)2 


17.65(5) 


225(6) 


16.5(4) 


17.4 


BaCu 2 Sb 2 


2.59(3) 


160(2) 


17.9(2) 


18.9 



Then the T dependence Pbg{T) / p RG (Q R ) is a function 
/ of the ratio y = T/G R and is given by Eqs. (21)-(23) 
as 



Ah 

/(y) = 4.226 259 y 5 / 
Jo 



(e x - 1)(1 - e~ x ) 



— dx. (24) 



Now we can fit the experimental p(T) data by the T 
dependence of the BG prediction according to 39 



p(T)=p Q +p(Q R )f(T/e R ), 



(25) 



where po, p(Qr) and Or are all independently adjustable 
fitting parameters. To carry out the fit of the experi- 
mental data by Eq. (25), we use an analytic Pade ap- 
proximant function (a ratio of two polynomials in y -1 ) 
that accurately fits f(y) in Eq. (24). The numerical 
coefficients in the approximant are obtained by a high- 
accuracy nonlinear-least-squares fit of numerical values 
of f(y) obtained from Eq. (24) by the Pade approximant 
as described in Ref. 39. 

The solid red curve in Fig. 8 is the least-squares fit 
of the p(T) data in Fig. 8 by Eqs. (24) and (25) using 
the Pade approximant function 39 in place of Eq. (24). 
The fit parameters obtained are po = 9.90(7) pflcm, 
p(& R ) = 34.9(6) /xO cm, and R = 245(4) K. The value 
of R is identical within the errors to the value ©d = 
246(1) K obtained from the fit of heat capacity data over 
the whole measured T range by the Pade approximant 
for the Debye function (Table IV). The value of TZ(& R ) 
can be calculated from the value of p(&r) using Eq. (23), 
which gives 7?.(0 R ) = 36.9 pilcm. The parameters ob- 
tained from the fit of p(T) for SrCu2As2 by the BG model 
are summarized in Table VII. 

To compare our experimental value of p(300 K) 
for SrCu2As2 with the value predicted by the Bloch- 
Griineisen theory, it is necessary to slightly modify the 
expression for 1Z(Q R ) in Eq. (22), which applies to 
monatomic metals, to allow application to polyatomic 
metals. Thus one obtains 39 



p BG (T = 6 R ) = 0.9 464 635ft(6 R ). 



(23) 



K 



7? 



7T 3 (37T 2 ) 1/3 ft 2 (J_ 

4n^afc B e R 



(26) 



12 



where the variables have the same meaning as before, 
except that 




is the average inverse mass of the atoms in the compound, 
n is the number of atoms per f.u. and Mi is the atomic 
weight of element i. In this case, the equivalent lattice 
parameter of a primitive cubic unit cell containing one 
atom is a = [V cc ii/ (nZ)} 1 / 3 (nZ = 10 here) because there 
are Z — 2 formula units of AC112X2, with n = 5 atoms 
per formula unit, per body-centered tetragonal unit cell 
with volume V ce n given in Table II. 

The number of current carriers per atom n ce ii in 
Eq. (26) can be estimated from the band theory value 29 of 
vp. Since m*vp = hkp, and using the one-band relation 37 
kp = (37r 2 n*) 1 / 3 where n* is the carrier density per band, 
one obtains 

n * = _L ( H£^L] 3 = 8.78 x 10 20 cm" 3 , (28) 
6-k z \ n J 

where to obtain the second equality we used the value 
of (vp) given in Eq. (9) and the value of m* given in 
Eq. (11). This gives 

n ce ii = A band n* % = 0.0327. (29) 
nZ 

Inserting this value of n ce n and the value 8r = 245 K 
from Table IV into Eq. (26) we obtain K(Q R ) = 
10.0 /zSl cm. The value of p(300 K) predicted by the 
BG model is then numerically calculated from Eq. (21) 
to be p(300 K) = 21.7 /ificm. The measured value 
p(300 K) = 54 ^ficm from Fig. 8 is a factor of 2.5 larger 
than this theoretical value. However, this should be con- 
sidered to be satisfactory agreement in view of the even 
larger factors of 3 to 8 discrepancy between experiment 
and theory for other sp-metals such as Au, Be, Cd, Zn, 
Mg, Al and In. 45 Such discrepancies can be attributed 
at least in part to deviations of the shapes of the Fermi 
surfaces from the spherical shape assumed in the Bloch- 
Griineisen theory. 

We conclude that both the magnitude and T depen- 
dence of the in-plane electrical resistivity of SrCu2 As2 are 
adequately fitted over the entire T range 1.8-300 K by 
the Bloch-Gruneiscn model in which the T-dependence 
of the resistivity arises solely from clectron-phonon scat- 
tering. 

V. PHYSICAL PROPERTIES OF SrCu 2 Sb 2 
CRYSTALS 

A. Heat Capacity 

The C P (T) data for SrCu2Sb2 are presented in Fig. 9 
for 1.8 < T < 300 K. The C P (T) data do not show any 




T(K) 

FIG. 9: (Color online) Heat capacity C p of a SrCu2St>2 single 
crystal versus temperature T. The solid red curve is a fit of the 
data by the sum of the contributions from the Debye lattice 
heat capacity CV Debye (T) and predetermined electronic heat 
capacity jT according to Eq. (12). Inset: Cp/T vs. T 2 below 
5 K. The straight line is the fit by C p /T = 7 + (3T 2 [Eq. (1)] 
for 1.8 K < T < 3.5 K. 



evidence for any phase transitions over this T range. Like 
SrCu 2 As 2 , the value of C p {T = 300 K) « 124 J/molK for 
SrCu2Sb2 is very close to the expected classical Dulong- 
Pctit value C v = 124.7 J/molK. A fit of the low- 
T C P (T) data below 3.5 K by Eq. (1) (solid line in 
the inset of Fig. 9) yields 7 = 3.9(2) mJ/molK 2 and 
/3 = 0.76(3) mJ/molK 4 . Using Eq. (5), the 7 value 
gives T>(Ep) = 1.65 states/eVf.u. for both spin direc- 
tions which is larger than that of SrCri2As2. The Debye 
temperature determined using Eq. (3) is Od = 234(3) K. 

A fit of the C P (T) data over the whole T range 1.8- 
300 K by the sum of an electronic contribution jT, with 
7 fixed at the above value, and the Debye lattice heat 
capacity contribution via Eqs. (12) and (13) gives 6d 
= 214(1) K which is somewhat smaller than the above 
value Op, = 234(3) K obtained from a fit to the low-T 
C P (T) data only. The solid curve in Fig. 9 is the fit of the 
data by the sum of electronic 7T and Debye lattice heat 
capacity terms where for the latter we used our Pade ap- 
proximant function. 39 The parameters obtained from our 
fits to C P (T) for SrCu2Sb2 are summarized in Table IV. 



B. Magnetization and Magnetic Susceptibility 

The ZFC anisotropic x(T) data of a SrCu2Sb2 single 
crystal in H = 3.0 T are presented in Fig. 10 together 
with the intrinsic \ values after correction for PM and 
FM impurity contributions as described below. Similar 
to the case of SrCu2As2, the \ has negative values. How- 
ever, the x anisotropy in SrCu2Sb2 is different from that 
of SrCu2As2- Here it is clear that Xc > Xab instead of 
Xc < Xab- The reason for this large qualitative difference 
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FIG. 10: (Color online) Zero-field-cooled magnetic suscepti- 
bility x = M/H of a SrCu2Sfc>2 single crystal versus temper- 
ature T in the temperature range 1.8-350 K measured in a 
magnetic field of 3.0 T applied along the c-axis (xc, H || c) 
and in the ab-plane (xab, H _L c) (solid symbols). The open 
symbols are the intrinsic susceptibility of SrCu2Sb2 after cor- 
recting for ferromagnetic and paramagnetic impurity contri- 
butions. 



in the magnetic anisotropy between the two compounds 
is not known in detail but is probably associated with 
the differences between the crystal structures of the two 
compounds. 

The anisotropic isothermal M(H) data for SrCu2Sb2 
at different T are presented in Fig. 11. The M(H) curves 
exhibit anisotropic diamagnetic behaviors with M C (H) > 
M ab (H), as expected from the \{T) = M(T)/H data 
in Fig. 10. The contributions from FM and PM impu- 
rities that cause the slight nonlinearities in the M (H) 
isotherms at low fields were estimated by fitting the 
M(H) data by Eqs. (14) and (15). The M s values were 
determined by fitting the M(H) data above 2.0 T and at 
T > 25 K by Eq. (14). These values were found to be in- 
dependent of T. The fitting parameters for both M a i,(H) 
and M C (H) isotherms at 1.8 K are listed in Table V where 
we found that good fits are obtained using 9 lmp = for 
both data sets. The fits by Eq. (15) are shown as the 
solid curves in Fig. 11. Similar to the case of SrCu2As2, 
both M s and the product fi mp Si nip have different values 
for H _L c and H || c, again indicating that the magne- 
tizations of the FM and PM impurities in SrCu2Sb2 are 
anisotropic. The M(H) isotherms at the other (higher) 



T were fitted using the values of /i mp and <Si m p obtained 
from the 1.8 K M(H) fits in order to extract the intrinsic 
anisotropic susceptibilities. We have corrected the mea- 
sured susceptibility data for these FM and PM impurity 
contributions, and the resulting intrinsic susceptibilities 
are shown by open symbols in Fig. 10. There are no 
significant differences between the uncorrected and cor- 
rected data above 50 K. For T < 50 K, the PM impurity 
Curie- Weiss-like contribution is nearly eliminated in the 
corrected data, yielding an intrinsic susceptibility that is 
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FIG. 11: (Color online) Isothermal magnetization M of a 
SrCu2Sb2 single crystal versus magnetic field H measured at 
the indicated temperatures for H applied (a) in the afe-plane 
(M ab , Hlc) and (b) along the c-axis (M c , H || c). The 
solid curves represent fits by Eq. (15) with the parameters in 
Table V, where (9i mp = was used for the fits. 



nearly independent of T, as expected. There is a 10% 
discrepancy between the measured value of Xab (Fig. 10) 
and that obtained from the high-field (H > 2 T) slopes of 
the M(H) isotherms for H _L c [Fig. 11(b)]. The reason 
for this discrepancy is not known, but it has no significant 
influence on our discussions or conclusions. 

The powder- and temperature-average of the intrinsic 
magnetic susceptibility is (%o) = —1-56 x 10~ 4 cm 3 /mol. 
The different contributions to this averaged intrinsic 
susceptibility were estimated in the same way as for 
SrCu2As2 in Sec. IV B above. The diamagnetic core sus- 
ceptibility is estimated as Xcore = —2.08 x 10~ 4 cm 3 /mol 
using the atomic diamagnetic susceptibilities. 42 From 
Eq. (19), the Pauli susceptibility is xp — 5.33 x 
10~ 5 cm 3 /mol using T>(Ep) — 1.65 states/eVf.u. for 
both spin directions determined above from the 7 value. 
The xl = —1-78 x 10~ 5 cm 3 /mol was obtained by tak- 
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FIG. 12: (Color online) In-plane electrical resistivity p of 
a SrCu2Sb2 single crystal versus temperature T. The solid 
curve is a fit by the Bloch-Griineisen model. 



ing m* — m c in Eq. (20). The Van Vleck suscepti- 
bility is then calculated from Eq. (18) to be xvv = 
1.66 x 10~ 5 cm 3 /mol. These magnetic susceptibility con- 
tributions are summarized in Table VI. 

The powder-averaged x of SrCu2Sb2 is significantly 
more negative than that of SrCu2As2, probably due at 
least in part to the larger Xcoro of Sb compared to that of 
As. On the other hand, the magnitude of the anisotropy 
of x is much larger in SrCri2Sb2 than in SrCu2As2, for 
reasons that are not clear to us. 



C. Electrical Resistivity 

The in-plane p(T) data for SrCu2Sb2 are presented in 
Fig. 12. The T dependence of p indicates metallic be- 
havior with a residual resistivity p — 13.6 pflcm at 
1.8 K and RRR « 3. The p(T) of SrCu 2 Sb 2 was ana- 
lyzed using the Bloch-Griineisen model. A least-squares 
fit by Eqs. (24) and (25) using our Pade approximant 
function 39 is shown as the red solid curve in Fig. 12. 
The fit parameters are po — 13.64(1) pQ, cm, /?(0r) = 
13.8(1) (tzOcm and 6 R = 153(1) K. Using Eq. (23), 
one obtains 7^(0r) = 14.6 pVi cm. These parameters 
are summarized in Table VII. 



VI. PHYSICAL PROPERTIES OF ORIENTED 

SrCu 2 (Aso.84Sb .i6)2 CRYSTALS 

A. Heat Capacity 

The C P (T) data for SrCu2(Aso.84Sbo.i6)2 are pre- 
sented in Fig. 13. Like SrCu2As2 and SrCu2Sb2, the 
C P (T) data of SrCu2(Aso.84Sb .i6)2 reveal no evidence 
of any phase transitions. The C P (T = 300 K) 
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FIG. 13: (Color online) Heat capacity C p of 
SrCu2(Aso.84Sbo.i6)2 versus temperature T. The solid 
red curve is a fit of the data by the sum of the contributions 
from the Debye lattice heat capacity CV Debye (r) and prede- 
termined electronic heat capacity 7T according to Eq. (12). 
Inset: C p /T vs. T 2 below 5 K. The straight line is a fit by 
Cp/T = 7 + /3T 2 [Eq. (1)] for 2.0 K < T < 3.5 K. 



122.8 J/molK is again close to the Dulong-Petit value 
CV = 124.7 J/molK. From a linear fit of C p /T vs. T 2 
data below 3.5 K by Eq. (1) (solid red line in the inset 
of Fig. 13), we obtained 7 = 2.3(2) mJ/molK 2 and /3 
= 0.83(2) mJ/molK 4 . The value of 7 yields V(E F ) = 
0.97 states/eVf.u. for both spin directions from Eq. (5) 
and the value of /3 gives 6 D = 227(2) K from Eq. (3). A 
value Gd = 237(1) K is obtained from a fit of the C P (T) 
data by the Debye lattice heat capacity model [Eqs. (12) 
and (13)] over the T range 1.8 < T < 300 K as shown 
by the solid red curve in Fig. 13. This value is similar to 
the value obtained from the low-T fit. The above specific 
heat parameters are compared in Table IV with those of 
the other three compounds studied here. 



B. Magnetization and Magnetic Susceptibility 

The anisotropic ZFC x( T ) data of 
SrCu2(Aso.84Sbo.i6)2 in H = 3.0 T are presented 
in Fig. 14. Like SrCu 2 As 2 and SrCu 2 Sb 2 , the x of 
SrCu2(Aso.84Sbo.i6)2 is diamagnetic. As can be seen 
from Fig. 14, Xab > Xo similar to the case of SrCu2As2 
but in contrast to that of SrCu2Sb2. Thus the anisotropy 
seems to be a characteristic of the crystal structure, 
since the structure of SrCu2(Aso.84Sbo.i6)2 is the same 
as that of SrCu2As2 but different from that of SrCri2Sb2. 

The anisotropic isothermal M(H) data for 
SrCu2(Aso.84Sbo.i6)2 at different T are presented 
in Fig. 15. The M{H) curves exhibit anisotropic behav- 
ior with M ab > M c , consistent with the x(T) = M(T)/H 
data discussed above. Furthermore, we observe that the 
M is almost proportional to H at each T, indicating 
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FIG. 14: (Color online) Zero-field-cooled magnetic suscepti- 
bility x = M/H of SrCu2(Aso.s4Sbo.i6)2 versus temperature 
T in a magnetic field H = 3.0 T applied along the c-axis 
(Xc, H || c) and in the ab-plane (Xab, H _L c). The solid stars 
represent the intrinsic susceptibility of SrCu2(Aso.84Sbo.i6)2 
obtained from the high-field slopes of M(H) isotherms. 



only very small amounts of FM and PM impurities in 
the sample. The intrinsic xo values obtained from the 
high-field (H > 2 T) slopes of the M(H) isotherms 
are shown as solid stars in Fig. 14. The powder- and 
temperature-average (xo) of Xo and t ne different contri- 
butions Xcorej XPj Xl and xvv to (xo) estimated in the 
same ways as in section IV are presented in Table VI. 
Interestingly, the (xo) of SrCu2(Aso.84Sb .i6)2 is not 
intermediate between those of SrCu2As2 and SrCu2Sb2. 



C. Electrical Resistivity 

The in-plane p(T) data for SrCu2(Aso.84Sbo.i6)2 are 
presented in Fig. 16. Like SrCu2As2 and SrCu2Sb 2 , the 
p{T) of SrCu2(Aso.84Sbo.i6)2 exhibits metallic behavior 
with po = 17.6 pQcm, which is the largest value among 
the four compounds studied here (see Table VII). The 
RRR « 2 is smaller than those of the end-point com- 
pounds SrCu 2 As 2 (RRR « 6) and SrCu 2 Sb 2 (RRR » 3) 
as expected from the atomic disorder on the As/Sb sub- 
lattices, and also of BaCu2Sb 2 (RRR « 15, see be- 
low) with yet another structure. The large value of 
po and the low RRR value for SrCu2(Aso.84Sbo.i6)2 are 
consistent with this compound being an sp-metal be- 
cause the conduction carriers evidently strongly sense 
the atomic disorder on the As/Sb sublattices as they 
hop along these sublattices. An analysis of the p(T) 
data of SrCu2(Aso.84Sbo.i6)2 using the Bloch-Grimeisen 
model in Sec. IV C (solid curve in Fig. 16) gives p = 
17.65(5) (itQcm, p(6 R ) = 16.5(4) fiQcm, 9 R = 225(6) K, 
and using Eq. (23), 1Z(Or) = 17.4 /zficm. These param- 
eters are compared in Table VII with those obtained for 
the other three compounds studied here. 
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FIG. 15: (Color online) Isothermal magnetization M of 
SrCu2(Aso.84Sbo.i6)2 versus magnetic field H measured at 
the indicated temperatures for H applied (a) in the afe-plane 
(M a i,, H _L c) and (b) along the c-axis (M c , H || c). 



VII. PHYSICAL PROPERTIES OF BaCu 2 Sb 2 
CRYSTALS 



A. Heat Capacity 

The C P (T) data for a BaCu2Sb2 crystal are presented 
in Fig. 17. The data reveal no evidence for any phase 
transitions. The data exhibit a saturation value of 
~ 125 J/molK near room temperature which is close 
to the expected classical Dulong-Petit value of Cy = 
124.7 J/molK. The origin of the weak upturn in C P (T) 
above 250 K is unknown and may be spurious, per- 
haps arising from a small inaccuracy in the high-T ad- 
denda calibration. We obtained the coefficients 7 = 
3.5(2) mJ/molK 2 and (3 = 0.65(2) mJ/molK 4 from a 
linear fit of the low-T C p /T vs. T 2 data below 3.5 K by 
Eq. (1) which is shown as the solid straight line in the 
inset of Fig. 17. The 7 value and Eq. (5) yield V(E F ) = 
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FIG. 16: (Color online) In-plane electrical resistivity p of 
SrCu2(Aso.84Sb>o.i6)2 versus temperature T. The solid curve 
is a fit by the Bloch-Griineisen model. 
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FIG. 17: (Color online) Heat capacity C p of a BaCu2St>2 sin- 
gle crystal versus temperature T. The solid red curve is a 
fit of the data by the sum of the contributions from the De- 
bye lattice heat capacity CvDobye(r) and predetermined elec- 
tronic heat capacity 7T according to Eq. (12). Inset: C p /T 
vs. T 2 below 5 K. The straight line is a fit by Eq. (1) for 
2.0 K < T < 3.5 K. 



1.49 states/eVf.u. for both spin directions. The f5 and 
Eq. (3) yield 6 D = 246(3) K. From a fit of the C P (T) 
data over the full T range by a sum of electronic and 
lattice contributions via Eqs. (12) and (13), we obtained 
a good fit using 0d = 204(1) K as shown by the solid 
red curve in Fig. 17. These specific heat parameters are 
summarized in Table IV. 



o 
E 

CO 

E 

o 

un 

o 



-4 
-6 
-8 
-10 
-12 



— 1 — i — | — 1 — i — 1 — 1 — | — 1 — i — i — r 

BaCu 2 Sb 2 
H 1 c, M/H 
H ± c, Intrinsic 
H || c, M/H 
H || c, Intrinsic 



_1_ 



50 100 150 200 250 300 350 
T(K) 



FIG. 18: (Color online) Zero-field-cooled magnetic suscepti- 
bility x = M/H of a BaCu2Sb2 single crystal versus tem- 
perature T measured in H = 3.0 T applied along the c-axis 
(Xc, H || c) and in the afe-plane (xab, H _L c) (solid sym- 
bols). The intrinsic susceptibilities after correcting for the 
small ferromagnetic impurity contributions are shown as open 
symbols. 



B. Magnetization and Magnetic Susceptibility 

The anisotropic x(T) data of a BaCu2Sb2 single crystal 
in H = 3.0 T are plotted in Fig. 18. Like the other 
three compounds, the susceptibility is diamagnetic over 
the whole T range. The data are weakly T-dependent. 
A weak anisotropy is seen below ^ 170 K with Xc > Xab 
as in SrCu2Sb2- 

The M (H) isotherms at different T are shown for 
BaCu2Sb 2 in Fig. 19. Consistent with the x(T) re- 
sults in Fig. 18, the M(H) curves of BaCu2Sb2 exhibit 
weakly anisotropic diamagnetic behavior with M C (H) > 
M a b(H)). The isotherms provide evidence for the pres- 
ence of FM impurities but no evidence for saturable PM 
impurities in the sample. In order to obtain the intrinsic 
X the M(H) isotherms at each T were fitted by Eq. (14) 
at high fields H > 2 T and the resulting T-independent 
FM impurity contributions (M s ofc = 0.15 Gcm 3 /mol and 
A/ S c = 0.20 Gcm 3 /mol) were subtracted from the M(T) 
measured at H = 3 T. The different values of M° b 
and M s c suggest an anisotropic FM impurity contribu- 
tion. To give the scale of the FM impurity magnetiza- 
tion, the above value M s ab = 0.15 G cm /mol is equivalent 
to the magnetization contribution from 12 molar ppm of 
Fe metal impurities. The intrinsic anisotropic suscepti- 
bilities obtained after subtracting the FM impurity con- 
tributions are shown by open symbols in Fig. 18. The 
corrections for the FM impurities are seen to be rather 
small. 

We have estimated the different contributions to the 
intrinsic x °f BaCu2Sb2 using the same methods as 
for the previous three compounds. The powder- and 
temperature-average over the T range 10 to 100 K is 
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FIG. 19: (Color online) Isothermal magnetization M of a 
BaCu2Sb>2 single crystal versus magnetic field H measured at 
the indicated temperatures for H applied (a) in the afe-plane 
(Mab, Hlc) and (b) along the c-axis (M c , H || c). 



(Xo) 



-8.02 x 10 5 cm 3 /mol, and we infer Xc 



-2.31 x 10~ 4 cm 3 /mol, xp = 4.81 x 10~ 5 cm 3 /mol, 
and xl = —1.60 x 10~ 5 cm 3 /mol (using m* — m c in 
Eq. 20). We thus obtain xvv = 1.19 x 10~ 4 cm 3 /mol 
from Eq. (18). A comparison of xvv with those of the 
other compounds (Table VI) suggests that xvv is largest 
in BaCu2Sb2- 



C. Electrical Resistivity 

The p{T) data for BaCri2Sb2 are presented in Fig. 20. 
The p(T) exhibits metallic behavior with po = 2.5 pflcm 
at 1.8 K and RRR « 15, indicating that this crystal 
has the highest perfection of the four compounds studied 
here. An analysis of the p{T) data in terms of the Bloch- 
Griineisen model in Eq. (25) yielded a good fit with po = 
2.59(3) pflcm, p(6 R ) = 17.9(2) pflcm, 6 R = 160(2) K, 
and using Eq. (23), TZ(& R ) = 18.9 /iOcm. The fit is 
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FIG. 20: (Color online) In-plane electrical resistivity p of a 
BaCu2Sb2 single crystal versus temperature T. The solid red 
curve is a fit by the Bloch-Griineisen model. 



shown as the solid red curve in Fig. 20. The fit param- 
eters are summarized in Table VII for comparison with 
those of the other samples. 



VIII. DISCUSSION 

The electronic and magnetic properties of SrCri2As2, 
SrCu2Sb2, SrCu2(Aso.84Sbo.i6)2 and BaCu2Sb2 indicate 
that they are all nonmagnetic sp-band metals. The mag- 
netic susceptibilities show temperature-independent dia- 
magnetism and the heat capacities show small Sommer- 
feld electronic specific heat coefficients, both indicative 
of small densities of states at the Fermi energy. The elec- 
trical resistivities show no evidence for a T 2 term that 
might indicate significant d-electron scattering, and in- 
stead are all well-described by the Bloch-Griineisen the- 
ory for scattering of sp-band electrons by longitudinal 
acoustic lattice vibrations. 

For SrCu2As2, our results are consistent with the band 
structure calculations for ThCr2Si2-type SrCu2As2 and 
BaCu2As2- 29 According to these calculations, the Cu 3d 
bands are centered about 3 eV below Ep and are com- 
pletely occupied, with little Cu 3d character in the bands 
at Ep. The same behavior occurs in Cu metal, where the 
Cu 3d bands are centered 3.4 eV below Ep with a total 
3d-band half- width of 1.4 eV. 46 The band calculations 
thus indicate that SrCu2As2 and BaCu2As2 are sp-band 
metals like Cu metal. 29 

In a chemical picture the Cu atoms in these com- 
pounds have a 3d 10 electronic configuration with a filled 
d-shcll and a formal oxidation state of Cu +1 , where the 
Cu 4s 1 valence electrons are itinerant. Assuming Cu +1 , 
Sr +2 and Ba +2 oxidation states in (Sr, Ba)Cu2As2 yields 
an unusual oxidation state of As~ 2 for the As atoms, 
in contrast to the oxidation states Fe +2 and As -3 in 
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(Ca,Sr,Ba)Fe2As2. 11 This in turn suggests the presence of 
a covalent bond between the As atoms in adjacent layers 
that are directly above/below each other along the c-axis 
[see Fig. 4(a)], i.e., As~ 2 = [As-As] ~ 4 /2. Such covalent 
interlayer pnictogen bonding has been known since the 
1980's for transition metal phosphides with the ThCr2Si2 
structure. 47-49 The variable degree of P-P covalent bond- 
ing in those compounds is reflected in a variable interlayer 
P-P distance <ip_p and a resulting variable c/a ratio. 50 
Notably, first-order transitions with increasing pressure 
between non-bonding and partially or fully bonding in- 
tcrplane P-P linkages were found for SrRli2P2 (Refs. 51, 
52) and SrNi 2 P 2 . 53 

In the recent field of iron arsenide and chalco- 
genide high-T c superconductivity, 11 Kreyssig, Goldman 
et al. 54,55 discovered a strong decrease in the c-axis lat- 
tice parameter and unit cell volume of ThCr2Si2-type 
CaFe2As2 in a first-order transition with increasing pres- 
sure p and dubbed this phase a "collapsed tetragonal" 
(cT) phase, as opposed to the ambient p (uncollapsed) 
tetragonal (T) phase. The T to cT transition occurs with 
increasing p at 0.4 GPa at low T, increasing to 1.7 GPa at 
300 K. 55 This first-order T to cT phase transition can be 
viewed as resulting from the discontinuous creation with 
increasing p of covalent As-As interlayer bonds, 55 as oc- 
curs at ambient p or with increasing p in some ThCi - 2Si2- 
type phosphides as discussed above. The "cT" terminol- 
ogy has been widely adopted by the community working 
on iron arsenide superconductors. 

The work by Kreyssig, Goldman et al. 54,55 moti- 
vated many additional experimental and theoretical stud- 
ies of the occurrence and properties of the cT phase 
in ThCr2Si2-type iron arsenides. A synchrotron x-ray 
diffraction study of SrFe2As2 at room T revealed a 
p-induced second-order transition at 10 GPa from the 
T to a cT structure for which electronic structure calcu- 
lations indicated the formation of As-As interlayer bond- 
ing as responsible for the lattice collapse. 56 Another syn- 
chrotron x-ray diffraction study at room T showed that 
a cT transition occurs in BaFe2As2 at an even higher 
p ~ 27 GPa, much higher than the transition p of 1.7 GPa 
found for the Ca compound at this T. A previous high- 
p synchrotron x-ray diffraction study of the same com- 
pound found that the onset p of the cT phase at room T is 
17 GPa under nonhydrostatic conditions and 22 GPa un- 
der hydrostatic conditions. 58 Synchrotron x-ray powder 
diffraction measurements versus T and p were also car- 
ried out for Cao.67Sr .33Fe2As2, with results intermediate 
between those obtained previously for the pure Ca and 
Sr end-members. 59 In the Ca(Fei_ x Co :I ;)2As2 system, the 
onset p for the cT transition at low T for x = 0.032, 0.051 
and 0.063 decreases with increasing a;. 60 At ambient p, 
a T-induced transition from the high-T T phase to the 
low-T cT phase occurs below ~ 90 K for x — 0.05-0.09 in 
the system CaFe2(Asi_ x P K )2. 61 The doped compounds 
Ca(Fei_2;Rh x )2As2 similarly exhibit an ambient pressure 
cT phase up to 50 K for x = 0.024 and up to 300 K for 
x = 0.19. 62 The T-p phase diagrams of lanthanide-doped 




FIG. 21: (Color online) Interlayer As-As distance c(as-As (a) 
and interlayer P-P distance dp_p (b) (left ordinates) and 
c/a ratio (right ordinates) versus 3d transition metal M for 
AM2A.S2 and AM2P2 compounds, respectively. The com- 
pounds with magnetically ordered ground states are enclosed 
in the green boxes. The compounds for which the magnetic 
ground states are not known are marked with "?". The ref- 
erences are given in Table VIII. 



Cai_ 3; i? 3 ;Fe2As2 (R = La, Ce, Pr, Nd) compounds, in- 
cluding T and cT phase regions, were also determined. 63 
In the following Sec. VIII A we examine more quantita- 
tively the systematics of interlayer P-P and As-As bond- 
ing in 3d transition metal phospides and arsenides with 
the ThCr2Si2 structure, and we will find that SrCii2As2 
is in the cT phase. Then in Sec. VIII B we discuss the 
relationship of the magnetic character of the ThCr2Si2- 
type 3c? transition metal arsenides and phosphides with 
the degree of interlayer pnictogen bonding and discuss 
how the properties of SrCu2As2 fit into this scheme. In 
Sec. VIIIC we discuss doping scenarios when Fe is par- 
tially substituted by other transition metals, including 
Cu, in the (Ca, Sr, Ba)Fe2As2 compounds. 



A. Pnictogen Interlayer Bonding 

Several estimates of the (single) covalent bond dis- 
tances dx-x have been made for X = As and P. In 
1988, Von Schnering and Honle inferred from a study of 
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TABLE VIII: Magnetic and structural properties of AM 2 X 2 compounds (^4 = Ca, Sr, Ba; M = Cr, Mn, Fe, Co, Ni, Cu; X = 
P, As) with the ThCr2Si2 structure, dx-x = (1 — 2zx)c and dx-x — {%zx — §) c are the interlayer and intralayer distances 
between the pnictogen atoms, respectively. Separate references are given for the crystal data and magnetic properties. For the 
latter, the abbreviation AFM means antiferromagnetic and SDW means itinerant spin density wave. 



Compound 


a (A) 


c(A) 


c/a 


zx 


dx-x (A) 


dT-x (A) 


Ref. 


Magnetic properties 


CaFe 2 As 2 (T) 


3.872(9) 


11.730(2) 


3.03 


0.3665(9) 


3.13(2) 


2.73(2) 


[5] 


AFM SDW [5] 


CaFe 2 As 2 (cT) 


3.9792(1) 


10.6379(6) 


2.67 


0.3687(7) 


2.79(2) 


2.52(2) 


[54] 


Fe moment quenched [54,55] 


CaCo2As2 


3.9831(2) 


10.2732(6) 


2.58 


0.3664(2) 


2.745(4) 


2.392(4) 


[66] 


AFM [66,67] 


CaNi 2 As 2 


4.053(6) 


9.90(2) 


2.44 


0.370 


2.574 


2.376 


[30,33] 


Pauli paramagnetic [33] 


CaCu 2 As2 


4.129(1) 


10.251(1) 


2.48 


0.3799(2) 


2.462(4) 


2.663(4) 


[68] 




SrCr 2 As 2 


3.918(3) 


13.05(1) 


3.33 


0.367 


3.471 


3.054 


[30,33] 


AFM [33] 


SrFe 2 As 2 


3.9289(3) 


12.3172(12) 


3.13 


0.36035(5) 


3.440(2) 


2.718(2) 


[69] 


AFM SDW [41] 


SrCo 2 As 2 


3.935(7) 


11.83(2) 


3.01 


0.362 


3.265 


2.650 


[30,33] 


Paramagnetic [70] 


SrNi 2 As 2 


4.1374(8) 


10.188(4) 


2.46 


0.3634(1) 


2.783(2) 


2.311(2) 


[71] 


No magnetic order [71] 


SrCu 2 As 2 


4.2725(1) 


10.2000(3) 


2.39 


0.3789(1) 


2.470(2) 


2.630(2) 




Nonmagnetic [This work] 


BaCr 2 As 2 


3.963(3) 


13.60(1) 


3.43 


0.361 


3.781 


3.019 


[30,33] 


Itinerant AFM [72] 


BaMn 2 As 2 


4.1686(4) 


13.473(3) 


3.23 


0.3615(3) 


3.732(8) 


3.004(8) 


[73] 


AFM [73] 


BaFe 2 As 2 


3.9633(4) 


13.022(2) 


3.29 


0.35424(6) 


3.797(2) 


2.715(2) 


[74] 


AFM SDW [40] 


BaCo2As2 


3.958(5) 


12.67(2) 


3.20 


0.361 


3.522 


2.813 


[30,33] 


Renormalized paramagnet [75] 


BaNi 2 As 2 


4.112(4) 


11.54(2) 


2.81 


0.3476(3) 


3.517(7) 


2.253(7) 


[76] 


Pauli paramagnetic [33] 


BaCu2As2 


4.446(5) 


10.07(1) 


2.26 


0.374 


2.538 


2.497 


[30,33] 


Pauli paramagnetic [33] 


CaFe 2 P 2 


3.855(1) 


9.985(1) 


2.59 


0.3643(3) 


2.710(6) 


2.283(6) 


[77] 


Pauli paramagnetic? [78] 


CaCo 2 P 2 


3.858(1) 


9.593(1) 


2.49 


0.3721(4) 


2.454(8) 


2.343(8) 


[77] 


Itinerant AFM [48,49] 


CaNi 2 P 2 


3.916(1) 


9.363(1) 


2.39 


0.3774(5) 


2.296(9) 


2.386(9) 


[77] 


Pauli paramagnetic [79] 


CaCui.75P2 


4.014(1) 


9.627(1) 


2.40 


0.3831(4) 


2.251(8) 


2.563(8) 


[77] 




SrFc 2 P 2 


3.825(1) 


11.612(1) 


3.04 


0.3521(8) 


3.43(2) 


2.37(2) 


[77] 


Pauli paramagnetic [80] 


SrCo 2 P 2 


3.794(1) 


11.610(1) 


3.06 


0.3525(5) 


3.42(2) 


3.38(2) 


[77] 


Stoner enhanced paramagnet [81] 


SrNi 2 P 2 


3.948(1) 


10.677(3) 


2.70 


0.3539(3) 


3.120(7) 


2.219(7) 


[82] 


Paramagnetic [83] 


SrCui. 75 P2 


4.166(1) 


9.607(1) 


2.31 


0.3805(4) 


2.296(8) 


2.507(8) 


[77] 




BaMn 2 P 2 


4.037(1) 


13.061(1) 


3.24 


0.3570(3) 


3.735(8) 


2.795(8) 


[77] 


AFM [84] 


BaFe 2 P 2 


3.840(1) 


12.442(1) 


3.24 


0.3456(4) 


3.84(1) 


2.38(1) 


[77] 


Pauli paramagnetic? [78] 


BaCo 2 P 2 


3.7994(3) 


12.391(1) 


3.26 


0.3461(3) 


3.814(8) 


2.382(8) 


[85] 




BaNi 2 P 2 


3.947(1) 


11.820(1) 


2.99 


0.3431(3) 


3.709(7) 


2.201(7) 


[82] 


Pauli paramagnetic [86] 



polyphosphides that dp_p = 2.23(7) A. 64 A recent survey 
in 2008 of the covalent radii of the elements by Cordero 
et al. gave d P _ P = 2.14(6) A and g? As -As = 2.38(8) A. 65 
As will be seen, the difference between the two dp_p val- 
ues as well as the error bars on any of these three dx-x 
values are negligible compared with the variations in the 
interlayer dx-x values across the various AM 2 (As,P) 2 
series of compounds. Here we will use the dp^p and 
c^As-As values of Cordero et al. for consistency. 

We have compiled in Table VIII crystallographic data 
from our work on SrCu2As2 together with data from the 
literature on ThCr 2 Si2-type AM 2 X 2 compounds (A = 
Ca, Sr, Ba; M = Cr, Mn, Fe, Co, Ni, Cu) with X = P 
and As. 5 < 33 < 40 ' 41 ' 54 < 55 > 66 - 86 The d X -x and the c/a ratio 
for X = As and P are plotted versus M in Figs. 21(a) 
and 21(b), respectively. The dp_p and c?As-As covalent 
bond distances of Cordero et al. 65 are plotted as the hor- 
izontal dashed lines in the respective figures. The cT 
phases in each figure are characterized by c?As-As or rfp_p 
values that approach the respective interlayer covalent 
bond distance. 

From Fig. 21, the d X -x and c/a values for the 
various ThCr2Si2-type compounds are strongly corre- 
lated, as previously documented by Jeitschko et al. for 



phosphides. 50 The c?As-As and dp-p are plotted versus 
c/a in Figs. 22(a) and 22(b) for the respective series of 
AM 2 X 2 compounds. An approximate linear relationship 
is seen between dx-x and c/a for each class of pnictides, 
with the phosphides following the linear behavior more 
precisely. Linear least squares fits of the data by 

dx-x =A + B- (30) 
a 

yield the parameters 

A = -0.28 A, 5 = 1.19 A, (X = As) (31) 
A = -1.78 A, 5 = 1.73 A. (X = P) 

The fits are shown by the green lines in Figs. 22(a) 
and 22(b), respectively. A plot of c?p_p versus c/a show- 
ing an approximate linear relationship for a more diverse 
range of ThCr 2 Si2-type phosphides is shown in Fig. 1 of 
Ref. 50. A plot of c?As-As versus c/a for ambient pressure 
(Ca, Sr, Ba)Fe2As2 and high-pressure cT CaFe2As2 show- 
ing a linear relationship was given in Fig. SI of Ref. 87. 

For SrCu 2 As 2 , our value c^as-As = 2.470(2) A is close 
to the As-As covalent bond distance 2.38(8) A, which 
confirms that the formal oxidation states of Cu and As 
atoms in this compound are Cu +1 and As -2 as discussed 
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FIG. 22: (Color online) (a) Interlayer As-As distance g!as-As 
and (b) interlayer P-P distance dp_p versus c/a in AM2AS2 
and AM2P2 compounds, respectively. Linear least-squares 
fits of the data in (a) and (b) by Eq. (30) are plotted as the re- 
spective green lines with parameters given in Eqs. (31). Filled 
symbols designate magnetically ordered compounds whereas 
the open symbols refer to compounds that do not have a mag- 
netically ordered ground state. The compounds for which the 
magnetic ground states are not known are marked with "?". 
The references are given in Table VIII. 



above. The +1 oxidation state of Cu is consistent with 
our p(T), C P (T) and x(T) measurements of SrCu 2 As 2 
as also previously discussed. In Sec. IX we suggest that 
the As-As interlayer covalent bonding in SrCu2As2 may 
be driven by the high stability of the Cu +1 d 10 electronic 
configuration rather than the Cu +1 oxidation state being 
a result of the As-As interlayer bonding. 



B. Correlations between Pnictogen Interlayer 
Bonding and Magnetic Properties 

The (Ca,Sr,Ba)M2AT2 compounds that exhibit magnet- 
ically ordered ground states are indicated in Fig. 21 by a 
box around them and in Fig. 22 by filled symbols. One 
sees that with the exception of cT phases CaCo2P2 and 
CaCo2As2 which are antiferromagnets, 48 ' 49 ' 66 ' 67 mag- 



netic ordering occurs only for the T phase compounds. 

Reehuis et al. 49 observed for a more extended set of 
AC02P2 compounds with the ThCr 2 Si2 structure, includ- 
ing A = rare earths, that either long-range magnetic or- 
dering of the Co atoms did not occur or FM ordering 
occurred if there is little or no P-P interlayer covalent 
bonding, whereas the character of the magnetic ordering 
changes to AFM if P-P interlayer bonding occurs (these 
are cT compounds in our terminology). The latter infer- 
ence is consistent with our datum for CaCo2P2 in Figs. 21 
and 22. 

Kreyssig, Goldman, Pratt et al. found that whereas 
CaFe2As2 exhibits SDW ordering in the T phase 
as previously reported, 5 the pressure-induced transi- 
tion to the cT phase quenches the magnetic order- 
ing as well as the Fe magnetic moment and associ- 
ated AFM fluctuations. 54 ' 55,88 Danura et al. found in 
the Ca(Fei_ a; Rh a: )2As2 system that after suppressing the 
SDW with increasing x, superconductivity appeared at 
about 13 K at 1 w 0.02, but then the superconductiv- 
ity was quenched at higher x when the compound went 
into the cT phase. 62 Furthermore, they found that con- 
verting a superconducting T sample with x w 0.02 to the 
cT phase under pressure destroyed the superconductivity. 
The high-p results of Uhoya et al. on BaFe2As2 suggested 
that the p-induced superconductivity with T c = 34 K at 
p = 1 GPa is suppressed in the higher-p stability range 
of the cT phase. 58 These results taken together suggest 
that the superconductivity in these phases is mediated 
by spin fluctuations. 

The degree of X — X covalent bonding and hence the 
formal oxidation states of the M atoms in the above 
AM2AS2 compounds are thus directly correlated with 
their magnetic and superconducting properties. 

From first principles calculations, Yildirim proposed 
early on that the degree of As-As bonding in CaFe2As2 
and the magnitude of the Fe spin are inversely 
related. 89,90 He suggested that the application of pres- 
sure on CaFe2As2 causes a reduction in the Fe moment 
which weakens the Fe-As bonding, and in turn, leads to 
As-As bonding and hence the cT phase. This prediction 
of a p-driven transition from a magnetic ground state in 
the T phase to a nonmagnetic one in the cT phase was 
confirmed by subsequent first-principles calculations (see, 
e.g., Ref. 91 in which the effects of nonhydrostatic pres- 
sure are also calculated) . Yildirim further suggested that 
the Fe moment should occur at p = even in the param- 
agnetic phase, 89 ' 90 which was confirmed by the observa- 
tion of AFM spin fluctuations above the Neel tempera- 
ture Tn by inelastic neutron scattering measurements on 
single crystals of the (Ca, Sr, Ba)Fe2As2 compounds. 11,92 



C. Transition Metal M Substitutions in 

Ba(Fei_ :c M :c )2As2 

Studies of the effects of partially substituting the Fe 
atoms in the ^4Fe2As2 compounds by other transition 
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metals showed that the long-range SDW ordering and 
associated lattice distortions in the undoped semimetal- 
lic parent compounds had to be largely suppressed before 
high-T c superconductivity could occur. 11 This was docu- 
mented, for example, by Canfield and colleagues who de- 
termined the T-x phase diagrams of Ba(Fei_ :c M :c )2As2 
(M = Co, Ni, Cu, Pd, Rh) from single crystal studies. 12 
They showed that the superconducting transition tem- 
perature T c (x), but not the AFM ordering tempera- 
ture Tn(x), approximately overlapped if the x-axis is 
replaced by the average d-electron concentration of the 
transition metal, thus suggesting the importance of this 
parameter to the superconducting properties. 12 Inter- 
estingly, they found that whereas Co and Ni substitu- 
tions for Fe induced bulk superconductivity, Cu substi- 
tutions did not even though they suppressed the struc- 
tural and SDW transitions as Co and Ni substitutions 
do. 12 This has been a conundrum in the Fe As-based 
superconductivity field. Angle-resolved photoemission 
spectroscopy (ARPES) measurements on some of these 
Ba(Fei_ x M 2; )2As2 systems support the contention that 
such aliovalent M-doping increases the electron concen- 
tration (see, e.g., the references in Ref. 93). 

On the other hand, as discussed above, applying pres- 
sure to the parent compounds can cause the same changes 
without doping. Furthermore, isovalent substitutions of 
Ru for Fe or P for As in BaFe2As2 suppress the SDW 
and induce high-T c superconductivity and thus have T-x 
phase diagrams similar to those of the Ba(Fei_ x M x )2As2 
compounds where M is an aliovalent transition metal. 11 
Remarkably, Dhaka et al. 93 found from angle-resolved 
photoemission spectroscopy (ARPES) measurements on 
Ba(Fei_2,Ru x )2As2 crystals that the electronic structure 
near E-p does not change with Ru doping. They specu- 
lated that the suppression of the SDW with Ru substi- 
tution and concomitant onset of superconductivity arise 
from a reduction of the Stoner enhancement factor of the 
static susceptibility at the ordering wave vector. 

Various theoretical studies have been carried out to de- 
termine whether the carrier concentration changes upon 
aliovalent transition metal substitutions for Fe in the 
AFe2As2 compounds. Density functional theory (DFT) 
calculations by Wadati et al. of the electronic structures 
of Ba(Fei_ x M x )2As2 compounds with M = Co, Ni, Cu, 
Zn, Ru, Rh and Pd suggested that these dopants do not 
change the carrier concentration because any additional 
ci-electrons are localized on the dopant sites. 94 On the 
other hand, another DFT study of disordered substitu- 
tions of 12.5% of Co for Fe in BaFc 2 As 2 by Berlijn et 
al. 95 found that a large chemical potential shift occurs 
due to the extra d-electron of Co as observed in ARPES 
experiments, 96-98 and that the chemical disorder induces 
important carrier scattering effects. They also found that 
the electron concentration due to the doping tends to pile 
up around the Co dopant sites as found by Wadati et 
al., 94 but that these electrons are associated with band 
states and are not localized. 95 

On the experimental side, local measurements of the 



electronic environments around the Fe and other atoms 
in Ba(Fei_a;M x )2As2 have been carried out. Fe K 
near edge x-ray absorption structure (NEXAS) mea- 
surements of the electron density around the iron site 
in Ba(Fei_ a; Co a ;)2As2 indicate that the extra doped 
ci-electrons are not concentrated at the Fe sites. 99 Simi- 
larly, Merz et al. found from L 2 ,3 NEXAS measurements 
on the Fe, Co and As atoms in the Sr(Fei_ x Co x )2As2 
system that no observable changes in the Fe, Co and As 
valences occur on Co-doping. 100 Resonant and nonres- 
onant x-ray absorption and emission spectroscopy mea- 
surements as well as core and valence level x-ray pho- 
toelectron spectroscopy measurements were carried out 
by McLeod et al. 101 on the Fe, Co, Ni and Cu atoms 
in Ba(Fei_; C M r )2As2 crystals with M = Co, Ni and Cu. 
They found that there is little charge transfer from the 
dopant atoms to the Fe atoms and that the iron arsenides 
arc not strongly correlated electron systems. They con- 
cluded that the Co and Ni substitutions provide addi- 
tional conduction electrons, but that the Cu 3d elec- 
trons are localized in the Cu 3d 10 shell and therefore that 
Cu-doping reduces the free carrier concentration at the 
Fermi level. 101 Resonant L3-edge Fe and Co photoclec- 
tron spectroscopy measurements of Ca(Fei^ a; Co a ;)2As2 
crystals with x — 0.056 by Levy et al. revealed that the 
Co 3d electrons participate in the formation of the Fermi 
surface, that the Fe and Co atoms have the same +2 oxi- 
dation state, and therefore that the extra d electron of Co 
is donated to the conduction band but whose charge den- 
sity is associated mainly with the Co dopant atoms. 102 
The totality of all these local measurements on Co and 
Fe indicate that the doped electrons of Co are itinerant 
and participate in the formation of the Fermi surface, 
but that their spatial densities are located near the Co 
atoms. 

Ideta et al. carried out a comprehensive ARPES 
study of Ba(Fei_ x M x ) 2 As 2 crystals with M = Ni (x = 
0.0375, 0.05, 0.08) and Cu (a; = 0.04, 0.06, 0.08) in- 
cluding measurements of both the in-plane afr-axis and 
out-of-plane c-axis Fermi surface dispersions and com- 
pared the results with previous ARPES measurements 
on Ba(Fei„ 2; Co 2 ;)2As2. 103 They found that volume en- 
closed by the electron Fermi surface(s) increased and the 
volume enclosed by the hole Fermi surface(s) decreased 
with increasing x for Co, Ni and Cu substitutions for Fe, 
as expected from electron doping in a rigid band picture. 
They also found that the (i-bands associated with the 
substituted Co, Ni and Cu atoms have increasing binding 
energy with respect to Ep , as expected from the increas- 
ing impurity potential. However, one also expects in a 
rigid band approach that the difference n e — between 
the numbers of electron and hole carriers per doped Co, 
Ni or Cu atom should be equal to the number of extra do- 
nated (i-electrons per atom. From the difference between 
the volumes enclosed by the electron and hole Fermi sur- 
faces, they found this prediction to be violated for Ni and 
Cu dopants, with the itinerant electrons doped per atom 
by Ni and especially by Cu smaller than the expected 
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values of two and three, respectively. 103 They suggested 
that n c — rij, is too small for Ni and Cu substitutions 
either because some of the donated electrons occupy im- 
purity bands that do not cross the Fermi level or they 
occupy localized states. When the authors plotted the 
SDW Neel temperature Tn and the superconducting T c 
versus n c — ?ih, they found universal behaviors for both 
quantities for all three dopants Co, Ni and Cu. Thus 
although the rigid band picture is not accurate with re- 
spect to electron doping, Ideta et al. have identified a 
quantity that controls both Tn and T c in electron-doped 
BaFe 2 As 2 . 103 



The Divergent Properties of SrCu2 As2 and Cu-doped 
BaFe2 As2 and the Role of As-As Interlayer Bonding 

In this paper, we have conclusively demonstrated from 
physical property measurements that SrCu2As2 is an sp- 
band metal. We have also demonstrated that As-As in- 
terlayer covalent bonding occurs in this compound, which 
is consistent with the Cu +1 oxidation and with the non- 
magnetic d 10 electronic configuration. In a band pic- 
ture, the Cu 3<i-bands at ~ 3 eV binding energy below 
E-p are fully occupied and the Cu e?-orbitals contribute 
little to the bands at Ep. 29 Therefore, on the basis of 
our results, the most natural and straightforward ex- 
pectation for the effect of partially substituting Cu for 
Fe in the (Ca,Sr,Ba)Fe2As2 compounds is that such Cu 
substitutions result in strong hole-doping instead of the 
expected 26 strong electron doping. In this case, the hole- 
doping resulting from replacing Cu for Fe in the AFe2As2 
compounds would be qualitatively similar to hole-doping 
resulting from partial substitutions of Mn or Cr for Fe. 
This analogous doping behavior would empirically ex- 
plain why strongly suppressed and non-bulk supercon- 
ductivity with T c < 2 K was observed in Cu-substituted 
BaFe2As2, 26 since partially substituting Mn or Cr for Fe 
does not induce superconductivity. 21-25 The origin of the 
latter behavior, however, is not understood microscopi- 
cally at present, especially since hole-doping by partially 
substituting K for Ba in BaFe2As2 induces superconduc- 
tivity with T c up to 38 K. 1 

However, as discussed above, it has been demonstrated 
both experimentally and theoretically that low-level sub- 
stitutions of Cu for Fe in BaFe2 As2 result in electron dop- 
ing rather than hole doping. Furthermore, from Figs. 21 
and 22, interlayer As-As bonding does not occur in 
BaFe2As2 whereas it does in BaCu2As2- These differ- 
ences indicate that the As-As bonding character and as- 
sociated Cu doping character in Ba(Fei_ a: Cu; C )2As2 must 
change with increasing x, as possibly suggested from 
thermoelectric power and Hall effect measurements on 
single crystals of this system. 105 At small x, the Cu evi- 
dently replaces the 3d 6 Fe +2 atoms as a 3c? 9 Cu +2 electron 
dopant, but at high concentrations would dope as 3d 10 
Cu +1 , which would correspond to a crossover between 
electron and hole doping with increasing x. However, to 
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FIG. 23: (Color online) (a) Unit cell volume V^ e ii versus com- 
position x in Ba(Fei_ x Cua;)2As2 (upper sets of data) and for 
the two endpoints x — and x = 1 of Sr(Fei_ a; Cu a ,)2As2. The 
Ba(Fei_ x Cu a; )2As2 data with blue diamond symbols are from 
Ni et al. 26 and the green square at x — 0.9 is a preliminary 
result of the present authors. The remaining data are from 
the literature. 11 All compounds have the bet ThCr2Si2-type 
structure. The two dashed lines represent Vegard's law for 
the respective series of compounds, (b) The c/a ratio versus 
x for the Ba(Fei_ a; Cu ;!; )2As2 system. The data symbols have 
the same meaning as in (a). 



accommodate the change in the interlayer As-As bonding 
from nonbonding to full covalent bonding and consider- 
ing the very dissimilar chemistry expected for Cu +1 and 
Fe+ 2 with increasing x, phase-separation (a miscibility 
gap) might occur in the Ba(Fei_ a; Cu :r ) 2 As 2 system, as we 
have demonstrated to occur in the Ba(Fei_ 2; Mn 2 ;)2As 2 
system. 104 Our preliminary studies of powder samples of 
Ba(Fei_ a; Cu a: )2As2 indeed suggest the presence of a mis- 
cibility gap in this system for 0.7<a;<0.9. 66 

Detailed structure refinements versus x of 
Ba(Fei_ x Cu K )2As2 would shed light on the nature 
of the As-As bonding and the doping character of Cu, 
but such structure refinements have not been reported 
up to now. However, Ni et al. 26 reported the lattice 
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parameters a and c and the unit cell volume V ce \i = a 2 c 
of Ba(Fei_ K Cu a: )2As2 versus x up to x = 0.35; their 
results for V^n versus x are shown in Fig. 23(a). Also 
included is our preliminary datum for a polycrystalline 
sample with x = 0.9 that was quenched from 900 °C, 
together with literature data for BaFe2As2, BaCu2As2, 
SrFe2As2 and SrCu2As2- 11 These four compounds have 
the same tetragonal ThCr 2 Si2-type structure. Between 
the x = and x = 1 end points of the Ba and Sr series, 
one sees from Fig. 23(a) that the volume decreases in 
each series by about the same amount of w 5 A 3 . How- 
ever, the data of Ni et al. 26 for Ba(Fei_ x Cu 2: )2As2 show 
a qualitative deviation from Vegard's law, exhibiting 
a strong increase in V ce \\ with increasing x instead of 
the expected decrease. This anomalous behavior may 
indicate interesting changes in the As-As bonding, Cu 
oxidation state and doping behavior versus x. In view 
of Fig. 22(b), we have plotted c/a versus x for the 
Ba(Fei_ x Cu :E )2As2 system in Fig. 23(b). The trend 
of the data suggest that interlayer As-As bonds may 
be beginning to form by x ~ 0.35. Detailed structural 
studies are clearly called for to clarify these characteris- 
tics versus Cu substitution concentration in FeAs-based 
ThCr2Si2-type systems. 



IX. SUMMARY AND CONCLUSIONS 

We have presented experimental data on the crys- 
tallographic properties in Sec. Ill and physical prop- 
erties in Sees. IV- VII of single crystals of SrCu2As2, 
SrCu2Sb2 and BaCu2Sb2 and of aligned clusters of 
SrCu2(Aso.84Sb .i6)2 crystals that were synthesized by 
the self-flux growth method. Contrasting structures 
were observed for SrCu2As2 (ThCr2Si2-type), SrCu2Sb2 
(CaBe2Ge2-type) and BaCu2Sb2 (a coherent intergrowth 
of ThCr 2 Si2-type and CaBe2Ge2-type unit cells, likely 
with a monoclinic distortion). The C P (T) and p(T) 
measurements show metallic behaviors for all four com- 
pounds. The x(T) data exhibit nearly T-independent 
diamagnetic behaviors. The \ °f SrCu2As2 is found to 
be larger in the a6-plane than along the c-axis, whereas 
the x values of SrCu2Sb2 and BaCu2Sb2 are larger along 
the c-axis. These differences in anisotropy appear to 
arise from the differences in the crystal structures be- 
tween these compounds. No evidence was observed from 
any of our measurements between 1.8 and 350 K for any 
phase transitions in any of the four compounds. 

In Sees. IV-VII, we estimated the different contribu- 
tions to the intrinsic \ and analyzed the C P (T) and p(T) 
data within the frameworks of the Debye model of lattice 
heat capacity and the Bloch-Griineisen model of resistiv- 
ity, respectively. The T-dependence of p is explained by 
electron-phonon scattering, as expected for an sp-band 
metal. In particular, there is no clear evidence for a T 2 
temperature dependence that would have indicated a sig- 
nificant contribution of d orbitals to the density of states 
at the Fermi energy Ep. Furthermore, our data rule out 



the possibility that the Cu ions are present as Cu +2 lo- 
cal magnetic moments with a 3d 9 electronic configuration 
and spin S = 1/2. 

Thus we conclude from the electronic and mag- 
netic properties of these four compounds that they 
are all nonmagnetic sp-metals, as previously predicted 
by D. J. Singh from electronic structure calculations 
for SrCu2As2 and BaCu2As2 with the ThCr2Si2-type 
structure. 29 These predictions are that the 3d bands of 
Cu are narrow and are centered about 3 eV below Ep, so 
they are completely filled, and there is little admixture of 
the 3d states into those at Ep. In chemical language this 
corresponds to a formal oxidation state of Cu +1 with a 
nonmagnetic 3d 10 electronic configuration. As discussed 
in Sec. VIII A, for SrCu2As2 the Cu +1 oxidation state is 
consistent with the presence of the observed interlayer co- 
valcnt As-As bonding, which together with the expected 
oxidation state Sr+ 2 yields As~ 2 = [As - As]~ 4 /2. 

On the other hand, for SrCu2Sb2 with the CaBe2Ge2- 
type structure, no Sb-Sb interlayer bonding is possible 
[see Fig. 4(b)] and the compound still shows the same 
physical properties as SrCu2As2, again indicating a Cu +1 
formal oxidation state and 3d 10 electronic configuration. 
It appears that the concept of formal oxidation state loses 
its meaning for Sb in SrCu2Sb2 because one cannot de- 
rive the Cu +1 oxidation state from an assumed oxidation 
state Sb -2 associated with the presence of Sb-Sb bond- 
ing. In both SrCu2As2 and SrCu 2 Sb 2 , the nonmagnetic 
character of Cu and the d 10 electronic configuration ev- 
idently arise because of the stability of this electronic 
configuration. In a band picture, the large Cu impurity 
potential causes the Cu 3d bands to be far below Ep, 
inhibiting their participation in the electronic conduc- 
tion. The above dichotomy between the As and Sb com- 
pounds in fact suggests that the As-As covalent bonding 
observed in SrCu2As 2 is driven by the stability of the 
Cu +1 oxidation state. A band structure calculation for 
SrCu2Sb2 and a comparison with that of SrCu2As2, to- 
gether with measurements of the local electronic configu- 
rations of Sr, Cu and Sb in both SrCu2As2 and SrCu 2 Sb 2 , 
would help to clarify these issues. 

As discussed in Sec. VIII C, a topic that is cur- 
rently being strongly debated is by how much, and even 
whether, transition metal substitutions on the Fe sites 
in j4Fe2As 2 compounds change the conduction carrier 
concentration. Early in the development of the iron ar- 
senide high-T c field, it appeared that substituting Co, 
Ni and Cu for Fe donates an additional one, two or 
three electrons to the conduction bands of the system, 
respectively, 12 but that expectation has been contested 
on several fronts. Our results on SrCu2As2 show that 
this compound is a collapsed tetragonal phase with Cu 
in the Cu +1 oxidation state with a d 10 electronic config- 
uration. This electron configuration suggests that sub- 
stituting Cu for Fe in (Ca,Sr,Ba)(Fei_ x Cu x )2As2 should 
result in hole-doping rather than the observed electron 
doping. Thus the electronic character of the Cu dopant 
and the strength of the As-As interlayer bonding are 
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both expected to drastically change between weakly Cu- 
substituted BaFc 2 As 2 and pure BaCu2As2, perhaps via 
a first-order lattice instability such as a miscibility gap 
in the Ba(Fe2- x Cu x )2As2 system. Detailed studies of 
the crystallography of (Ca,Sr,Ba)(Fci_. E Cu x )2As2 versus 
x would help to clarify the changes in the As-As inter- 
layer bonding and the nature of the carrier doping by Cu 
with increasing x. 
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